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ABSTRACT 


An  experimental  investigation  to  measure  particle  size  distributions  in  the 
plume  of  sub-scale  solid  rocket  motors  was  conducted.  A  phase-Doppler  particle 
analyzer  (PDPA)  in  conjunction  with  three-wavelength  extinction  measurements 
were  used  in  a  specially  designed  particle  collection  probe  in  an  attempt  to 
determine  the  entire  plume  particle  size  distribution.  In  addition,  a  laser  ensemble 
particle  sizer  was  used  for  comparative  data.  The  PDPA  and  Malvern  distributions 
agreed  in  the  observed  modes  near  1  and  4.5  |im  diameter  {d).  Scanning  Electron 
Microscope  (SEM)  picmres  of  collected  particles  were  in  good  agreement  with  the 
measured  Malvern  Sauter  Mean  diameter  of  2.59  ^im.  Data  analysis  indicate 
that  less  than  3%  of  the  total  mass  of  the  particles  was  contained  in  particles  with 
diameter  d  <  0.5  pm.  Therefore,  the  PDPA,  which  can  typically  measure  particles 
down  to  a  minimum  diameter  of  0.5  pm  with  a  dynamic  range  of  50:1, 

can  be  used  by  itself  to  determine  the  particle  size  distribution.  Multiple 
wavelength  measurements  were  found  to  be  very  sensitive  to  inaccuracies  in  the 
measured  transmittances. 
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I.  INTRODUCTION 


Solid  rocket  motors  have  used  aluminum  as  a  metallic  fuel  in  composite  and 
composite-modified  double-base  propellants  for  many  years.  Typically,  5  -  20  %  of  the 
mass  of  the  propellant  is  contained  in  crystalline  aluminum  particles  mixed  with 
approximately  70  %  ammonium  perchlorate  (AP)  as  an  oxidizer.  The  remaining  percent 
of  the  mass  is  contained  in  binders,  plasticizers,  or  burning  rate  catalysts.  Metallic  fuels 
help  performance  as  measured  by  an  increased  specific  impulse  7^^.  When  metallic  fuels 
combust,  they  react  with  the  oxidizer  and  form  metal  oxides.  For  aluminum,  the  product 
oxide  is  AI2O3,  also  called  alumina.  These  compounds  are  initially  in  the  liquid  phase  at 
the  high  temperatures  in  the  motor  chamber,  but  then  become  solid  with  the  temperature 
drop  through  the  nozzle  and  into  the  plume. 

In  aluminized  propellants  the  physical  process  of  combustion  has  been  under 
investigation  for  over  thirty  years.  However,  the  generally  accepted  theory  is  as  follows: 
An  AP/binder  flame  starts  the  process.  The  20  -  30  |xm  aluminum  particles  melt  at  933  K 
on  the  propellant  surface  and  sinter  together  forming  much  larger  agglomerates  but  an 
oxide  shell  forms  around  the  agglomerate  and  inhibits  further  combustion  until  the 
oxide’s  melting  point  is  reached  at  2327  K.  These  agglomerates  detach  the  surface  and 
enter  the  flame  zone.  There  the  oxide  shell  melts,  allowing  inner  aluminum  to  vaporize 
as  well  as  form  more  oxide.  An  oxide  "cap"  forms  on  one  side.  Thus,  clouds  of  burning 
agglomerates,  burning  unagglomerated  droplets,  and  product  oxide  droplets  exist  in  the 
motor.  When  the  droplets  bum,  the  vaporized  aluminum  forms  smoke  oxide  (diameters  < 
2  |xm)  that  streams  around  the  agglomerate.  The  larger  droplets  with  their  oxide  cap  form 
5-150  |xm  residual  oxide.  Thus,  any  particle  size  distribution  in  the  motor  is  expected  to 
bebimodal.  [Ref.  1-3] 

As  the  particles  pass  through  the  nozzle  various  processes  occur.  Collision 
coalescence  can  cause  larger  particles  to  form,  whereas,  shattering  produces  just  the 
opposite.  It  has  been  both  measured  and  modeled  that  the  larger  particles  entering  the 
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exhaust  nozzle  cannot  make  the  rapid  turn  at  the  throat,  resulting  in  their  being 
concentrated  more  along  the  plume  centerline.  The  smaller  particles  are  capable  of 
turning  with  the  gas  flow  and  have  been  observed  throughout  the  plume  cross-section. 
Losses  occur  due  to  two-phase  flow  of  liquid  or  solid  particles  and  gas  in  the  nozzle. 
These  can  be  further  defined  in  terms  of  thermal  and  velocity  lag  losses  which  vary  with 
the  size  of  particles.  Thus,  knowledge  of  the  particle  size  distribution  is  important  to 
determine  these  losses. 

Besides  determining  losses,  there  are  other  reasons  why  particle  size  information  is 
desired.  Plume  signature  and  flowfield  codes,  such  as  the  SIRRM  and  SPF,  require  this 
information.  Tactically,  the  AljOj  in  a  missile  plume  makes  the  plume  very  opaque  and 
visible,  therefore  leading  to  missile  detection,  tracking,  and  potentially  targeting.  For 
example,  when  viewing  a  Space  Shuttle  launch  from  several  miles  away,  the  plume 
emission  of  the  Shuttle  Rocket  Booster  (SRB)'s  is  highly  visible  but  the  huge  columns  of 
aluminum  oxide  exhaust  are  unforgettable.  Size  distribution  information  is  also 
important  for  determining  the  available  particle  damping  in  the  motor,  slag  formation, 
heat  transfer,  and  nozzle  flow  computations  [Ref.  1].  For  example,  Blomshield  [Ref.  4] 
conducted  dry  combustion  bomb  experiments  with  SRB  propellants  containing  different 
manufacturer’s  AP  to  evaluate  the  variation  of  the  resulting  AljOj  size  distribution 
between  them  and  their  impact  on  slag  formation  in  the  SRB's.  This  work  was  done  as  a 
result  of  the  pressure  anomaly  during  STS-54. 

The  properties  of  aluminum  oxide  are  important  for  use  with  different  particle 
sizing  techniques  and  calculations.  The  optical  properties  of  alumina,  such  as  the  index 
of  refraction,  are  important  for  use  in  optical  techniques  as  well  as  the  previously 
mentioned  computer  codes.  In  the  plume,  both  liquid  and  solid  alumina  can  exist  due  to 
the  range  of  temperature.  The  density  of  AljOj  varies  with  temperature  and  phase.  Reed 
[Ref.  2,  p.9]  lists  the  following  two  formulae  for  liquid  and  solid  alumina  respectively: 

Pl  =  2.98 -(1.13x10-3  XT’- 2327)  liquid  (1.1) 
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ps=4.00-(1.16xl0-^)7’  solid 


(1.2) 


where  p  is  density  in  g/cc  and  T  is  temperature  in  K.  Salita  [Ref.  3,  Figure  13]  provides  a 
plot  of  the  variation  of  aluminum  and  alumina  densities  with  temperature  as  well  as  their 
ratio.  Regarding  the  complex  index  of  refraction  (m  =  n-  ik),  Reed  [Ref.  2]  reviewed 
others'  results  indicating  a  range  for  the  AljOj  n  from  1.65  to  1.76.  The  complex  part  k, 
which  is  the  index  of  absorption,  has  shown  high  variability  but  is  almost  zero  for  the 
solid  phase  in  the  plume.  The  liquid  phase  has  a  much  higher  k  around  0.5  - 1.0  x  lO"^  . 

The  Naval  Postgraduate  School  (NPS)  has  been  conducting  research  in  the  area  of 
solid  rocket  motor  particle  sizing  both  in  the  motor  and  in  the  plume  using  various 
techniques  for  several  years.  Optical  techniques  have  received  the  most  attention  since 
they  offer  the  potential  for  obtaining  particle  information  without  disturbing  the  flow,  i.e., 
they  are  non-intrusive.  However,  collection  probes  are  sometimes  required  to  limit  the 
sample  volume  because  the  plume  is  too  opaque  for  forward  scattering  measurements. 
Various  optical  techniques  have  been  applied  individually  in  different  regions  of  the 
motor  and  plume.  Each  technique  has  a  limited  particle  size  measurement  range.  The 
present  investigation  used  a  phase-Doppler  particle  analyzer  (PDPA)  in  conjunction  with 
a  three-wavelength  extinction  measurement  in  an  attempt  to  determine  the  entire  plume 
particle  size  distribution.  In  addition,  a  laser  ensemble  particle  sizer  was  used  for 
comparative  data.  The  phase-Doppler  and  three-wavelength  extinction  measurements 
were  made  simultaneously  in  a  particle  collection  probe.  A  composite  propellant  with 
20%  aluminum  was  used  in  a  sub-scale  motor  to  provide  the  particle-laden  plume. 
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II.  BACKGROUND 


A.  PARTICLE  SIZE  DISTRIBUTIONS 

Particle  size  is  usually  represented  by  the  diameter  (in  microns)  of  an  equivalent 
size  spherical  particle.  For  a  collection  of  particles,  size  information  is  often  given  in  the 
form  of  a  size  or  frequency  distribution  where  frequency  is  plotted  against  particle 
diameter. 

There  are  various  types  of  frequency  distributions  depending  upon  the  parameter 
of  interest.  Typically,  this  is  number  or  mass  and  the  distribution  is  termed  either  a 
number  or  mass  distribution.  Frequency  is  another  loose  term  that  can  represent  several 
things.  For  example,  a  number  distribution  could  use  frequency  to  mean  actual  number, 
normalized  number,  or  number  percent  for  a  given  diameter. 

Due  to  experimental  limitations,  frequency  is  determined  empirically  for  diameter 
intervals  where  diameter  D  is  a  continuous  random  variable.  Other  terms  synonymous 
with  diameter  intervals  are  classes,  class  intervals,  or  size  bins.  The  middle  diameter  in 
each  class  is  sometimes  called  the  class  mark. 

Classification  of  data  leads  to  representation  using  a  special  type  of  bar  chart 
known  as  a  histogram.  It  uses  adjacent  rectangles  where  the  base  of  each  rectangle  is 
along  the  abscissa  and  represents  the  minimum  to  maximum  diameters  of  the  particles  in 
that  class.  The  height  of  each  rectangle  corresponds  to  the  ordinate  and  represents  the 
frequency  of  the  particles  in  that  class.  For  histograms  the  frequency  is  usually  given  per 
diameter.  For  example,  a  number  distribution  would  have  as  frequency  the  number 
percent  in  each  class  divided  by  the  class  width.  The  area  of  each  rectangle  then 
represents  the  number  percent  in  the  class.  This  is  a  very  useful  form  because  a 
probability  density  function  (PDF)  p(d)  can  be  fit  to  the  histogram. 

Synonymous  terms  for  probability  density  function  are  probability  function, 
density  function,  or  distribution  function.  The  PDF  can  be  visualized  as  the  smooth  curve 
resulting  from  connecting  the  mid-points  of  the  tops  of  the  rectangles  as  the  widths  of  the 


5 


rectangles  (class  widths)  approach  zero.  The  area  under  a  PDF  equals  1  or  100%, 
meaning  the  probability  P  that  the  diameter  is  between  -«>  and  +<>=  is  100%.  Symbolically 
P(-«>  <D<  +oo)  =  1,  or  mathematically 

f2p{D)dD=l  (2.1) 

Since  there  are  no  particles  with  diameters  less  than  zero,  the  lower  limit  can  be 
set  to  zero.  Similarly,  the  probability  P  that  the  diameter  is  between  two  arbitrary  sizes  a 
and  b  is  just  the  area  under  the  curve  from  a  to  Z?  or 

?{a<D<b)  =^\%{D)dD  (2.2) 

When  the  lower  limit  of  the  integral  above  is  zero  then  the  integral  represents  the 
probability  P  that  the  diameter  is  less  than  b.  This  has  a  special  name;  the  cumulative 
probability  function  or  cumulative  distribution  function/.  For  a  distribution  of  particles, 
replacing  b  with  d,  this  function /(^f)  represents  the  percent  (or  fraction)  of  particles  less 
than  diameter  d.  Mathematically, 

m  =  P(D  <  cO  =  p{D)  dD  (2.3) 

B.  NUMBER,  SURFACE/AREA,  AND  VOLUME/MASS  DISTRIBUTIONS 

Some  particle  sizing  techniques  count  particles  and  therefore  determine  a  number 
distribution  directly.  Others,  such  as  sieving,  separate  particles  based  upon  their  surface 
area,  volume,  or  weight.  Thus  it  is  paramount  to  know  what  distribution  a  given 
technique  will  yield  and  to  be  able  to  compare  these  different  distributions.  Additionally, 
the  total  number,  surface  area,  volume,  mass,  or  weight  of  the  particles  in  the  distribution 
are  used  for  comparison  and  analysis.  When  actual  values  are  of  interest  as  opposed  to 
percent  (e.g.,  number  instead  of  number  %),  then  the  distribution  function  often  differs 
from  the  PDF  by  a  factor  of  the  total  parameter  of  interest  (e.g.,  total  number  of  particles 
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in  the  distribution).  Take,  for  example,  the  number  distribution  N(d)  =  N\  PDF  where  N 
is  the  total  number  of  particles  in  the  distribution.  The  PDF  is  the  number  percent  per 
diameter,  and  N(d)  is  the  number  per  diameter.  Integration  of  N(d)  over  all  D  results  in 
N.  Integration  of  N(d)  from  to  yields  the  number  of  particles  between  and  ^2- 
Since  the  surface  area  of  a  sphere  is  it  cF,  the  surface  or  area  distribution  S(d) 
corresponding  to  N(d)  is  S(d)  =  N(d)  xnS  where  S{d)  is  the  area  per  diameter. 

Similarly,  the  volume  distribution  V{d)  corresponding  to  N{d)  is  V{d)  =  N{d)  xitd^l  6 
where  V(d)  is  the  volume  per  diameter.  For  a  constant  particle  density,  the  mass 
distribution  looks  the  same  as  the  volume  distribution  and  differs  numerically  by  a  factor 
representing  the  density.  Corresponding  area  and  volume  or  mass  PDF's  can  be  found  by 
dividing  S(d)  and  V(d)  by  S  and  V  the  total  area  and  volume  of  particles  in  the 
distribution,  respectively.  The  constants  S  and  V  can  be  found  by  integrating  S(d)  and 
Vid)  over  all  D.  Since  S(d)  and  V(d)  only  differ  from  their  corresponding  PDF's  by 
constants  they  have  the  same  shape  as  the  PDF's,  just  different  values. 

C.  S AUTER  MEAN  DIAMETER 

Various  diameters  are  used  to  compare  distributions.  They  are  referred  to  as 
mean,  average,  or  equivalent  diameters.  Just  as  the  different  distributions  are  equivalent 
to  N(d)  times  d  raised  to  some  power,  the  mean  diameters  can  be  represented  for  a 
continuous  distribution  as 


dpq  - 


j°^N(D)DPdD 


or  for  a  discrete  distribution  as 

I,N{d)dPM 
~  j:N{d)d^M 


(2.4) 


(2.5) 
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where  N(d)  is  the  number  per  class  width  Ad  [Ref.  5],  One  widely  used  mean  diameter  is 
the  Sauter  or  Volume-Surface  mean  diameter 

D.  LOG-NORMAL  DISTRIBUTION  FUNCTION 

In  certain  cases,  such  as  ground  particles,  the  histogram  has  a  Gaussian  or  normal 
distribution  shape  when  the  frequency  is  plotted  versus  the  logarithm  of  the  diameter. 
When  plotted  on  a  linear  scale  the  distribution  is  skewed  or  stretched  to  the  right.  Such  a 
distribution  is  said  to  be  log-normal. 

Like  a  Gaussian  or  normal  distribution,  the  log-normal  distribution  can  be 
characterized  by  two  parameters,  the  first  quantifying  the  central  tendency  and  the  second 
determining  the  spread  about  this  center.  For  a  Gaussian  distribution,  these  two 
parameters  are  the  mean  diameter  d^  (or  arithmetic  mean  J,,,)  and  the  standard  deviation 
c.  Similarly,  for  a  log-normal  distribution  these  parameters  are  the  geometric  mean 
diameter  d^  and  the  geometric  standard  deviation  a^.  Two  other  diameters  are  sometimes 
used  to  describe  central  tendency;  the  median  and  the  mode.  The  median  splits  the 
distribution  in  two  by  number  of  samples.  In  probability  terms,  if  all  the  N  sampled 
particles  are  put  in  order  of  increasing  size  and  assigned  a  number,  called  the  rank  (from 
one  to  AO,  then  the  median  diameter  is  that  of  the  particle  with  rank  (N  +  l)/2.  On  the 
other  hand  the  diameter  with  the  highest  frequency,  i.e.,  the  peak  of  the  distribution,  is 
termed  the  mode  diameter.  The  number  of  peaks  can  be  specified  by  referring  to  the 
distribution  as  monomodal,  bimodal,  or  trimodal.  It  is  important  to  note  that  d^  is  also 
the  median  and  the  mode  diameter  for  a  normal  distribution.  However,  for  a  log-normal 
distribution  d^  is  also  the  median  diameter,  but  is  neither  the  mode  nor  the  mean.  For  a 
log-normal  distribution: 


dM  =  dgexp(-ln  Gg) 

(2.6) 

dm  =  i/gexp(0.5  In^Gg) 

(2.7) 
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where  is  the  mode  diameter  and  is  the  mean  diameter  [Ref.  6]. 

For  a  log-normal  distribution  of  particles,  the  other  representative  distributions 
(i.e.,  number,  surface  area,  or  volume/mass)  are  also  log-normal  with  the  same 
However,  each  distribution  has  a  different  d^.  To  differentiate  between  them  it  helps  to 
use  an  additional  subscript.  For  example,  d^^  refers  to  a  number  distribution  whereas  d^^ 
refers  to  a  mass  distribution.  In  some  sources,  count  with  a  subscript  c  is  used  instead  of 
number.  The  Hatch-Choate  Equations  relate  the  various  d^'s.  One  of  these  is 

dgn  —  dgftt^^l^  ^  3  In  (cT^) J  (2.8) 


The  following  relation  is  also  helpful  in  relating  the  various  mean  diameters  to  d^. 


dpq  —  dgm^^]p 


(2.9) 


— —  ^/^^exp 


In  {Cg) 


(2.10) 


Since  different  distributions  are  of  interest  to  different  researchers,  rather  than  discussing 
several  d^  's  it  is  easier  to  use  a  mean  diameter  such  as  djj.  Note:  one  aspect  of  d^2 
counterintuitive.  As  one  might  expect  for  a  given  dg„,  as  increases  the  number 
distribution  is  skewed  more  to  the  right  and  J32  increases.  However,  if  djj  is  held  constant 
and  Gg  increased  d^^  decreases  and  the  number  distribution  shifts  left,  becoming  less 
skewed  to  the  right.  The  mass  geometric  mean  d^^,  however,  still  increases  and  the 
volume/mass  distribution  is  more  skewed  to  the  right. 

Sometimes  log-normal  distributions  are  plotted  on  log-probability  paper  since y(i/) 
plots  as  a  straight  line  with  slope  a^.  The  geometric  mean  d^  is  readily  identified  since  it 
corresponds  to  an  ordinate  value  of  0.5  or  50%.  All  the  distributions,  i.e.,  number,  area, 
and  volume,  plot  as  parallel  lines.  If  the  distribution  is  bimodal  rather  than  monomodal, 
then  it  has  an  s-shape  between  the  two  lines  making  up  the  individual  modes.  The 
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drawback  to  the  log-probability  plot  is  its  unavailability  in  standard  plotting  software 
such  as  spreadsheets.  According  to  the  MathWorks  however,  such  a  plotting  function  is 
available  for  MATLAB  in  their  Statistics  Toolbox. 

Several  forms  of  the  log-normal  distribution  function  are  presented  in  the 
literature.  Three  forms  are  listed  below: 


pid)  = 


exp 


Inidldg) 

2' 

_  J^Cg  _ 

(2.11) 


p{d)  =  -^ - 

■yjln  logic g)d 


f 

logid/dg) 

2" 

logic g)_ 

(2.12) 


pid)  = 


1 

lniCg)d 


exp 


Inid/dg)  ' 

2' 

_  ln(G^)_ 

(2.13) 


Some  sources  use  a  version  of  one  of  the  three  equations  above  without  the  d  in  the 
denominator.  [Ref.  6-10] 

It  is  therefore  important  to  determine  the  correct  form.  Three  criteria  are  used  to 
check  them: 

1 .  Integration  of  the  PDF  (area  under  the  curve)  from  if  =  0  to  +«>  must  equal  one 
(i.e.,y(+oo)=l). 

2.  Usingy(if),  the  ratio  of  d^^  (d  where  fid)  =  0.84  or  84%)  over  d^  id  where 
fid)  =  0.50  or  50%)  must  equal  c^. 

3.  The  function  must  behave  as  depicted  in  references  for  different  values  of  d 

and  G^.  * 

MATLAB  can  be  used  for  analysis.  The  functions  QUAD  and  QUAD8  performed 

numerical  integration,  although  an  exact  solution  for  fid)  is  available  for  Equations  (2.1 1) 

and  (2.13)  as  described  in  the  next  section.  Figure  1  shows  a  plot  of  Equations  (2.1 1)  to 

(2.13).  All  three  equations  without  the  d  in  the  denominator  fail  criterion  1.  However,  it 

should  be  noted  that  some  sources  integrate  the  PDF  with  respect  to  ln(d)  vice  d.  This  is 
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the  same  as  integrating  a  PDF,  with  a  c/  in  the  denominator,  with  respect  to  d  since 
</[ln(D)]  =  {\ID)dD  in  differential  form.  As  is.  Equation  (2.12)  fails  criterion  1  as  is 
readily  apparent  by  inspecting  the  figure.  Equation  (2.1 1)  fails  criterion  2.  Equations 
(2.1 1)  and  (2.12)  fail  criterion  3.  Equation  (2.13)  passes  all  criteria  and  thus  was  used  in 
the  analytical  work.  Note  however  that  Equation  (2.1 1)  is  a  PDF  and  fits  in  as  Kerker 
[Ref  6,  p.  357]  describes  "a  family  of  logarithmically  skewed  distributions".  [Ref  6] 


Figure  1.  Plot  of  Equations  (2.1 1),  (2.12),  and  (2.13)  as  PDF  #1,  #2,  and  #3 
respectively. 

E.  EXACT  SOLUTION  FOR  THE  LOG-NORMAL  F(D) 

The  general  form  of  the  log-normal  distribution  function  is 


/jN  1  ln(ii)-a 

p{d)  =  -= —  exp  K - =4= — 

j^^d  yip 


(2.14) 


where  p  and  a  are  constants  [Ref.  10].  Most  often,  for  solid  rocket  motor  particle  sizing 
P  =  In(ag)  and  a  =  In(iZ^)  as  in  Equation  (2.13).  Although  numerical  integration 
techniques  can  be  used  to  obtain  an  approximate y(<i),  there  is  an  exact  solution  using  the 
error  function  erf(x). 


erf(x)  =  ^  Jq  exp  {-f)dt 
and  using  Equations  (2.3)  and  (2.14) 


T'  1 

ln(D)-al^ 

[■ 

.  yip  J  ’ 

dD 


By  letting 


t  = 


ln(D)-a 

V2P 


then 


dt  = 


1 


■dD 


yipD 

and  the  limits  of  integration  are 
ln(D)-a 


t=  lim 

D-^d 


L  V2P  . 


In(c0  -  g 

y2p 


upper  limit 


(2.15) 


(2.16) 


(2.17) 


(2.18) 


(2.19) 


t=  lim 

£>-40 


ln(£))-a 

V2P 


lower  limit 


(2.20) 


Therefore, 
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(2.21) 


ln(<f>-o 


Jit 


(2.22) 


(2.23) 


(2.24) 


(2.25) 


(2.26) 


F.  MULTIPLE-WAVELENGTH  EXTINCTION  TECHNIQUE 

Cashdollar  [Ref.  12]  discusses  the  use  of  light  extinction  measurements  for 
extracting  particle  size  distribution  parameters.  Light  entering  a  cloud  of  particles  is 
absorbed,  scattered,  and  transmitted.  Transmittance  T  is  the  ratio  of  transmitted  intensity 
to  incident  intensity.  Light  extinction  is  the  sum  of  absorption  and  scattering.  When  Mie 
theory  is  applied  to  light  scattered  by  single  spherical  particles  the  extinction  coefficient 
Q  can  be  determined  as  a  function  of  wavelength  X,  particle  diameter  d,  and  complex 
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index  of  refraction  m  —  n  -  ik(k  is  the  absorption  index).  Bouguer's  law  provides  the 
transmittance  through  a  monodisperse  (single  size)  distribution  of  particles: 


T  =  exp 


^3QCn,L^ 

^  J 


(2.27) 


where 

T  =  transmittance 

Q  =  dimensionless  extinction  coefficient 

~  mass  concentration  of  particles  (mass  of  particles  per  volume  of  gas) 
L  =  path  length  through  the  cloud  of  particles 

p  =  particle  density 

d  =  particle  diameter 


The  volume  concentration  of  particles  (volume  of  particles  per  volume  of  gas)  can  be 

related  to  C„  by 


To  obtain  accurate  values  for  extinction,  it  is  necessary  to  have  the  detector  not  see 
forward  scattered  light.  An  approximation  for  the  required  detector  viewing  half  angle  0^ 
to  satisfy  this  requirement  is  [Ref.  12] 


01  <7.0- 
2  d 


(2.29) 


Dobbins  revised  the  Bouguer  transmission  law  to  apply  to  a  poly  disperse  system  (many 
sizes)  [Ref  12]; 


2pdi2  J 


(2.30) 


where  Q  is  the  average  extinction  coefficient,  found  for  a  continuous  distribution  by 
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(2.31) 


J”  Q{P)N{D)D^dD 
\°°^N{D)D'^dD 


or  for  a  discrete  distribution  by 


i:Q{d)N{d)cFM 

'LN(d)<fM 


(2.32) 


Since  Q  depends  upon  the  distribution,  then  for  a  log-normal  distribution  it  is  a  function 
of  J32,  "k,  and  m.  Taking  the  natural  logarithm  of  both  sides  of  Equation  (2.30)  and 
ratioing  two  equations  for  different  wavelengths  results  in 


\nT(ki) 

lnT(k2) 


Qiki) 


(2.33) 


This  is  the  basis  of  the  multiple-wavelength  technique.  Transmittances  for  different 
wavelengths  of  light  passing  through  the  same  volume  of  particles  over  identical  path 
lengths  are  experimentally  measured.  A  Mie  code  is  used  to  generate  0's  for  the  set 
wavelengths  and  for  various  values  of  J32,  and  m.  Then,  using  Equation  (2.33)  the 

ratios  can  be  compared  and  further  iterations  performed  as  necessary  to  determine  the 
best  values  of  and  m.  The  more  wavelengths  used  the  more  ratios  to  compare  and, 
presumably,  the  more  accurate  the  method.  The  number  of  ratios  for  a  set  of  n 
wavelengths  is  equal  to  the  sum  of  integers  less  than  n  (e.g.,  for  six  wavelengths  there  are 
1+2  +  3+  4  +  5=15  ratios).  Note  that  the  number  of  independent  ratios  is  only  n  - 1 
(e.g.,  five  for  six  wavelengths).  The  accuracy  of  the  method  is  also  increased  by 
spreading  the  wavelengths  across  the  spectrum  [Ref.  13].  To  determine  d^^’ 

m  =  n-ik(4  variables)  a  minimum  of  four  independent  ratios  are  required.  Figures  2,  3, 
and  4  show  plots  of  the  Mie  code  0's,  0's,  and  0  ratios  for  the  three  wavelengths  used  in 

the  experiments  and  an  assumed  log-normal  distribution.  Notice  from  Figure  4  that  as  J32 
increases  the  0  ratios  flatten  out.  For  this  reason  the  technique  does  not  work  for 
distributions  with  lots  of  mass  in  large  particles. 
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QBAR 


QBAR  RATIO 


Figure  4.  Average  Extinction  Coefficient  Ratios  for  Three  Wavelengths  and  a 
Log-Normal  Distribution. 
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III.  EQUIPMENT  DESCRIPTION  &  SETUP 


A.  ROCKET  MOTOR 

For  all  of  the  test  firings,  a  sub-scale  motor  was  used  as  shown  in  Figure  5.  The 
motor  consists  of  a  stainless  steel  cylindrical  chamber  with  a  2"  inside  diameter,  0.625 
wall  thickness,  and  overall  length  of  9.25".  The  head  of  the  chamber  is  closed  off  by  a  1" 
thick  end  plate  with  rubber  o-ring  seal  and  secured  by  six  bolts.  The  nozzle  end  of  the 
chamber  has  a  slightly  larger  inside  diameter  to  allow  the  nozzle  insert  to  slide  in  place, 
and  is  sealed  by  another  end  plate  and  six  bolts.  A  simple  converging-diverging  nozzle 
insert  is  used.  Gas  leakage  around  the  nozzle  is  prevented  by  the  use  of  two  o-rings  fit 
into  grooves  around  the  exterior  of  the  nozzle  insert.  A  pressure  transducer  and  pressure 
relief  fitting  are  screwed  into  the  chamber.  Prior  to  motor  firing,  a  hydraulic  dead-weight 
tester  is  used  to  determine  the  calibration  slope  and  zero  for  the  pressure  transducer.  The 
pressure  relief  fitting  consists  of  a  steel  elbow  with  a  1000  psia  burst  disk.  The  propellant 
is  bonded  to  the  head-end  plate  and  inner  chamber  wall  using  Permatex  High  Temp  RTV. 
The  igniter  consists  of  a  3/8"  steel  bolt,  described  below,  that  is  screwed  into  the  optimum 


De  =  0.365" 


20%  Al  propellant 


Not  shown:  pressure  transducer  &  1000  psia  relief  valve 


! 

/ 


probe  tip  locations 


Figure  5.  Sub-Scale  Motor. 
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of  four  available  holes  angled  toward  the  head  of  the  chamber.  The  remaining  holes  are 
filled  with  standard  bolts.  All  o-rings  are  lubricated  with  Dow  High  Vacuum  Grease. 

1.  Propellants 

Two  different  propellants  provided  by  the  Air  Force  Phillips  Laboratory  were 
used.  A  1"  thick  and  approximately  2"  diameter  propellant  grain  was  cut  from  the  stock 
using  a  large  knife,  then  the  desired  cutouts  were  made  using  a  three-piece  steel  cutting 
tool  and  a  press.  The  tool  could  make  0.6"  I.D.,  1.4"  O.D.,  1.5"  I.D.  and  2"  O.D.  cutouts. 
The  grain  geometry  was  chosen  to  optimize  the  burning  surface  area  for  a  neutral  bum 
with  a  chamber  pressure  of  at  least  250  psia  lasting  for  about  one  second.  RTV  was  used 
to  bond  each  grain  in  the  motor  as  well  as  to  inhibit  certain  surfaces  from  burning. 

For  testing  of  the  experimental  setup,  especially  the  three-wavelength  apparatus, 
several  firings  were  conducted  using  a  propellant  (AC- 14)  that  was  on  hand.  Its 
composition  is  shown  in  Table  1.  It  was  undesirable  for  data  analysis  due  to  the  silicon 
which  caused  oxides  of  silicon  or  aluminum  together  with  mullite  in  the  exhaust.  To 
achieve  the  desired  performance  several  grain  configurations  were  tried.  Initially,  a  single 
1"  thick  end-burning  grain  (2"  O.D.)  was  tried,  then  two  1"  thick  grains  were  tried,  one 
end-burning  (1.4"  O.D.),  the  other  radial-  and  end-buming  (1.5"  I.D.,  2"  O.D.).  Finally, 
good  performance  was  obtained  using  two  1"  thick  radial-  and  end-  burning  grains.  The 
first  had  a  2"  O.D.  and  a  1.5"  I.D.  (0.25"  web).  RTV  was  used  to  inhibit  the  downstream 
end  and  bond  the  grain  to  the  chamber  wall  just  downstream  of  the  igniter  hole.  The 
second  grain  was  cut  from  the  center  of  the  first.  It  had  a  1.4"  O.D.  and  a  0.6"  I.D.  (0.4" 
web).  RTV  was  used  to  bond  the  upstream  surface  to  the  endplate  and  to  inhibit  the  outer 
radial  surface. 
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ID 

AC-14 

Composition:  AP 

67.15% 

A1 

12.0% 

R45M 

10.14% 

Si 

6.0% 

Dioctyl  Adipate 

3.91% 

IPDI 

0.78% 

Triphenyl  Bismuth 

0.02% 

Density  (Ibm/in^) 

Pp  =  0.06 

Burning  rate  pressure  coefficient 

a  =  0.0086 

Burning  rate  pressure  exponent 

n  =  0.566 

Characteristic  exhaust  velocity  (ft/s) 

c*  =  4957 

Table  1.  Test  Propellant  AC- 14  Specifications. 


As  previously  mentioned  aluminized  propellants  have  been  the  subject  of  particle 
sizing  for  the  past  thirty  years.  Therefore,  it  was  desired  to  obtain  such  a  propellant  since 
various  parameters  are  better  known  (e.g.,  particle  index  of  refraction)  and  this  enables 
comparison  to  past  research.  The  USAF  Phillips  Lab  provided  four  propellants,  all  of 
which  were  class  1.3,  conventional  composites.  The  one  selected  for  this  thesis  research 
was  PS-1.  Its  composition  (by  mass  %)  and  characteristics  are  listed  in  Table  2.  Initially, 
a  single  1 "  thick  end-burning  grain  (2"  O.D.)  was  tried;  however,  the  required  small 
nozzle  throat  diameter  and  propellant's  high  aluminum  content  resulted  in  clogging  of  the 
throat.  A  single  perforated  end-  and  radial-burning  grain  was  then  chosen  and  used  for 
the  remainder  of  the  research.  It  was  1"  thick  with  a  2"  O.D.  and  a  0.625"  I.D. 
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ID 

PSl 

Composition:  AP  (coarse) 

50%  200|xm 

AP  (fine) 

20% 

A1 

20% 

HTPB 

10% 

Density  (Ibm/in^) 

pp  =  0.0667 

Burning  rate  pressure  coefficient 

a  =  0.0328 

Burning  rate  pressure  exponent 

n  =  0.376 

Characteristic  exhaust  velocity  (ft/s) 

c*  =  5119 

Table  2.  Operational  Propellant  PS- 1  Specifications. 


2.  Igniter 

Each  pyrotechnic  squib  igniter  was  processed  as  follows:  First,  a  hollow  3/8" 
steel  bolt  had  two  holes  drilled  in  the  bolt  head  to  accomodate  the  lead  wires.  Two 
approximately  8"  long  copper  lead  wires  were  then  cut  and  some  insulation  was  stripped 
off  both  ends  with  an  X-Acto  knife.  The  wires  were  passed  through  the  holes  and  then  a 
small  piece  of  nickel-chromium  wire  was  soldered  (using  ruby  flux  and  pure  solder)  to 
the  ends  leaving  about  a  1/16"  gap.  Next,  the  wire  bridge  was  pulled  inside  the  cavity  in 
the  bolt  and  the  bolt  head  was  epoxied  on  the  outside  to  secure  the  wires.  After  the  epoxy 
had  dried  the  cavity  was  filled  with  BKNO3  and  a  paper  hole  from  a  hole  punch  was 
glued  on  the  end  with  Elmer's  Glue.  Continuity  checks  were  performed  at  various  times 
in  the  process  to  verify  that  the  wires  were  not  shorted  by  touching  the  inside  of  the  bolt. 

For  firing  the  igniter,  lead  wires  were  connected  to  a  12  volt  battery  with  battery 
charger  and  a  firing  button  in  the  control  room.  When  the  firing  circuit  was  completed  by 
depressing  the  button,  the  heat  dissipated  in  the  nickel-chromium  bridge  wire  flashed  the 
BKNO3.  The  burning  matter  burst  through  the  paper  hole  and  impacted  the  grain,  thereby 
initiating  propellant  combustion. 
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3.  Nozzle 


A  copper  nozzle  insert  was  selected  to  obtain  reasonable  chamber  pressure.  The 
following  steady-state  mass  conservation  equation  was  used: 


’(Ah' 

\  ppClC  V\c* 

1  gc 

where 

=  chamber  pressure 
Aj  =  burning  surface  area 
A,  =  nozzle  throat  area 
Pp  =  propellant  density 
a  =  burning  rate  pressure  coefficient 
c*  =  characteristic  exhaust  velocity 

=  characteristic  exhaust  velocity  efficiency 
=  gravitational  constant  (for  English  units) 
n  =  burning  rate  pressure  exponent 

The  MICROPEP  computer  program  [Ref.  14]  was  run  with  the  propellant  ingredients  to 
obtain  and  c*.  See  the  output  in  Appendix  A.  Based  upon  previous  results,  the 
efficiency  T]^,  was  assumed  to  be  0.9.  Then,  the  available  nozzle  throat  areas  were 
substituted  to  calculate  p^.  Due  to  clogging  of  the  throat  from  the  high  aluminum  content, 
the  largest  A,  that  would  yield  p^  >  250  psia  was  selected.  The  copper  nozzles  were 
converging-diverging  and  Table  3  shows  the  selection  that  was  available. 

For  initial  testing  with  propellant  AC- 14  several  nozzles  were  used.  The 


d,  (in) 

0.207 

0.340 

0.239 

0.165 

0.452 

0.216 

0.270 

0.200 

0.208 

0.268 

0.286 

d,  (in) 

0.378 

0.440 

0.365 

0.230 

0.662 

0.395 

0.808 

0.368 

0.425 

0.500 

0.322 

e  =  AA 

3.335 

1.675 

2.332 

1.943 

2.145 

3.344 

8.956 

3.386 

4.175 

3.481 

1.268 

Con/Di 

(deg) 

45/15 

45/15 

45/15 

45/15 

45/15 

45/15 

??? 

??? 

??? 

45/15 

??? 

Table  3.  Nozzle  Selection  Where  d,  is  the  Nozzle  Throat  Diameter,  is  the  Nozzle 
Exit  Diameter,  e  is  the  Nozzle  Area  Ratio,  and  Con/Di  are  the  Nozzle  Converging  and 
Diverging  Half-Angles  Respectively. 
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d,  =  0.165"  nozzle  was  tried  with  the  single  end-buming  grain  and  resulted  in  twice  the 
expected  p^,  probably  due  to  clogging.  Next,  the  first  two-grain  configuration  was  paired 
with  the  d,  =  0.239"  nozzle  and  resulted  in  a  progressive  initial  bum  and  a  long  low 
tail-off.  Finally,  the  second  two-grain  configuration  was  paired  with  the  d,  =  0.268" 
nozzle  with  good  results. 

For  the  PS-1  firings,  two  nozzles  were  used.  The  d,  =  0.207"  nozzle  was  tried 
with  the  single  end-buming  grain  but  clogging  caused  a  highly  progressive  bum.  For  use 
with  the  single  end-  and  radial-burning  grain,  the  d,  =  0.239"  nozzle  was  chosen  and 
Figure  6  shows  the  expected  ideal  performance.  The  actual  performance  resulted  in  a 
higher  chamber  pressure  of  about  350  psia.  Post-firing  inspection  of  the  nozzle  revealed 
a  fairly  uniform  layer  of  AljOj  reducing  d^  to  about  0.15". 


Figure  6.  Expected  Ideal  Pressure-Time  Trace  for  PS-1 . 


B.  PROBE  AND  TUBES 

A  particle  collection  probe  that  consists  of  a  narrow,  hollow  steel  box  with  a 
conical  tip  for  swallowing  a  portion  of  the  plume  exhaust  was  previously  developed  [Ref. 
20].  The  design  of  the  probe  was  to  swallow  the  normal  shock,  forming  weaker,  internal 
oblique  shocks  that  would  not  break  up  the  particles.  It  was  positioned  to  avoid  Mach 
disk  locations.  Windows  on  either  side  had  been  used  with  optical  sizing  techniques  and 
a  filter  on  the  downstream  end  had  been  used  to  collect  particles  for  scaiming  electron 
microscope  (SEM)  analysis.  For  the  current  uses,  the  windows  were  removed  and  steel 
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tubes  added  on  the  sides  to  enclose  the  PDPA  beams  and  to  vertically  enclose  the  Hg-Xe 
beam  described  below.  The  insides  of  the  tubes  were  spray-painted  flat  black.  To  protect 
from  bumthrough,  an  additional  steel  plate  was  welded  on  the  funnel  side  tube.  To 
ensure  that  the  PDPA  lenses  would  not  be  coated  by  exhaust,  circular  plate  glass  was 
taped  in  the  ends  of  the  side  tubes.  A  nitrogen  line  with  a  sonic  choke  was  connected  to 
the  probe  and  provided  an  ejector  flow  around  the  probe  tip  where  it  exhausted  into  the 
probe  body.  The  back  pressure  in  the  probe  was  regulated  to  provide  and  approximate 
isokinetic  sampling  condition.  For  some  runs  particles  were  wiped  from  the  inside  of  the 
box  for  SEM  analysis.  RTV  was  put  on  parts  of  the  front  of  the  probe  body  for  thermal 
protection.  See  the  diagrams  in  Appendix  B. 

C.  WHITE  LIGHT  SOURCE  AND  DIODES 

An  Oriel  Model  66002  Hg-Xe  arc  lamp  with  a  f/1  condensing  lens  and  200  W 
power  supply  was  used  as  the  light  source.  It  was  mounted  below  the  lab  table  and  a 
magnetically  held  mirror  was  positioned  to  reflect  the  beam  up  through  a  1/2"  diameter 
hole  in  the  table.  From  the  top  of  the  table,  tubing  enclosed  the  beam  up  into  the  probe 
and  above  the  probe  to  the  diode  box.  Both  tubes  were  spray-painted  flat  black  inside. 
The  lower  tube  had  a  0.5"  O.D.  and  7/16"  I.D.  and  a  collar  adapted  it  to  the  probe's  lower 
vertical  tube  O.D.  A  larger  upper  tube  (1"  O.D.,  7/8"  I.D.)  was  used  to  cut  down  on 
internal  reflection  and  was  fit  over  the  probe's  upper  vertical  tube.  The  diode  box 
contained  an  adjustable  pinhole,  two  beamsplitters,  three  EG&G  Model  HUV-4000B 
diodes,  each  with  a  Melles  Griot  narrow  band  interference  filter  (0.01  pm  bandwidth)  to 
measure  wavelengths  0.3130, 0.5461,  or  0.8500  pm,  neutral  density  filters,  and  an 
NPS-built  preamplifier  for  the  diodes.  The  pinhole  and  neutral  density  filters  were 
adjusted  to  get  as  high  an  output  voltage  as  possible  but  less  than  10  volts  (LABTECH 
Notebook  limitation;  the  diodes  saturated  near  15  volts).  The  total  length  from  the 
condensing  lens  to  each  diode  was  about  80".  From  the  probe  center  to  the  diode  pinhole 
was  about  36".  Therefore,  since  the  pihole  was  open  to  about  d=0. 19"  the  detector 
viewing  half-angle  was  about  0.15°,  easily  satisfying  the  requirement  of  Equation  (2.29). 
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For  one  run  a  light  chopper  was  used  between  the  arc  lamp  and  mirror.  Figure  7  (lower 
curve)  shows  the  spectral  irradiance  from  0.200  to  0.900  pm  of  the  light  source  and  the 
bandpass  filter  locations. 

D.  PDPA  SYSTEM 

The  Aerometrics  Phase  Doppler  Particle  Analyzer  (PDPA)  system  consists  of  a 
laser  transmitter,  receiver,  doppler  signal  analyzer  (DSA),  and  computer  digital  signal 
processor  (DSP)  and  software.  A  Lexel  Model  95-4  argon  ion  laser,  tuned  to  0.5145  pm 
and  operating  at  2  W,  with  beam  waist  adjustor,  polarization  rotator,  and  steering  mirrors 
provided  the  input  beam  to  the  transmitter.  The  transmitter  uses  a  Bragg  cell  to  split  the 
incoming  laser  beam  into  two  beams  of  equal  intensity,  a  Zero  order  beam  and  a  40  MHz 
shifted  First  order  beam,  then  bends  the  beams  with  prisms  and  mirrors  and  converges 
them  with  the  exit  lens.  The  beams  converge  until  they  cross,  forming  an  ellipsoidal 
Gaussian  volume,  approximately  10"  in  front  of  the  transmitter.  When  a  particle  passes 
through  the  crossing  volume,  light  is  scattered.  For  particles  with  d  >  0.5  pm  this  is 
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scattering  in  the  Mie  and  geometric  optics  regimes  of  scattering.  To  simplify  processing 
the  system  assumes  only  the  geometric  optics  regime  applies.  The  receiver,  consisting 
essentially  of  optics  and  photo  multiplier  tubes,  is  positioned  at  a  predetermined  optimum 
angle  to  collect  a  portion  of  the  scattered  light.  Aerometrics  ran  simulations,  given 
estimated  particle  indices  of  refraction  for  different  motor  locations,  to  determine  this 
optimum  angle  for  a  linear  phase  response  over  the  range  of  particle  diameters.  They 
found  that  a  scattering  angle  of  130°  (backscatter  or  reflection)  was  best  for 
through-the-chamber  particle  sizing  and  a  scattering  angle  of  45°  (forward  scattering  or 
refraction)  was  best  for  through-the-plume  measurements  as  presented  in  this  paper.  The 
receiver  was  mounted  on  a  support  at  a  45°  angle  and  approximately  8"  from  the  crossing 
volume.  The  horizontal  tubes  of  the  probe  were  designed  to  fit  around  the  collars  of  the 
transmitter  and  receiver.  [Ref.  16] 

Data  acquisition  and  processing  was  done  with  a  personal  computer  (PC)  and  the 
DSA.  The  DSP  was  located  in  the  PC  and  was  used  to  quickly  perform  discrete  and  fast 
Fourier  transforms  on  the  data.  The  DSA  performed  several  functions.  It  provided  the 
input  to  the  transmitter  to  create  the  two  beams  and  conducted  calibration  and  data 
processing.  The  software  created  the  user  interface  through  menus  for  setup,  calibration, 
data  acquisition,  processing,  display,  and  storage.  The  PDPA  was  externally  triggered  to 
take  data  through  the  LABTECH  Notebook  program  run  on  a  separate  PC  as  described 
below.  [Ref.  16] 

The  PDPA  software  rejects  particles  for  various  reasons.  The  most  important 
ones  for  this  paper  are  multiple  particles  and  undersize  particles.  Only  single  particles 
can  be  accurately  accounted  for,  therefore  if  the  system  detects  multiple  particles  that 
sample  is  rejected.  Particles  with  diametersas  small  as  approximately  0.2  pm  are  detected 
but  if  less  than  0.5  pm  they  are  rejected  by  the  system.  No  information  is  available  for 
particles  less  than  0.2  pm  since  the  scattered  light  intensity  is  too  low  to  be  detected.  The 
PDPA  has  Probe  Volume  Correction  (PVC)  that  allows  for  correcting  the  measured 
distribution  for  the  smaller  particles.  However,  at  least  100  valid  samples  must  be  taken 
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Figure  8  shows  the  experimental  setup  for  the  PDF  A  and  three-  wavelength  extinction 
measurements  in  the  probe. 


Figure  8.  Equipment  Setup  for  PDPA  and  Three- Wavelength  Measurements  (Deflector, 
Video  Camera,  and  Supporting  Hardware  Not  Shown). 


E.  DEFLECTOR 

A  pneumatically  actuated,  computer  controlled  steel  deflector  plate  was  positioned 
between  the  motor  nozzle  and  probe.  Its  purpose  was  to  protect  the  probe  and  only  allow 
exhaust  entry  into  the  probe  during  steady  burning  (i.e.,  the  plateau  of  the  chamber 
pressure-time  trace  for  a  neutral  burning  grain).  Some  firings  burned  holes  through  the 
plate,  requiring  repair  by  welding  and  reinforcement. 

F.  DATA  ACQUISITION  SOFTWARE  &  HARDWARE 

A  486DX/33  MHz  PC  controlled  the  experiment  through  a  Windows™ 
application  known  as  LABTECH  Notebook.  Icon  blocks,  representing  analog  and  digital 
input  or  output  devices,  with  connecting  arrows  were  used  in  the  software  to  visualize  the 
control  of  equipment  and  the  flow  of  data.  See  Figure  9.  Each  icon  block  could  be 
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opened  to  set  various  parameters.  When  the  data  acquisition  file  was  run  the  following 
events  were  executed: 

1 .  The  diode  and  pressure  transducer  outputs  were  recorded  for  the  entire  run  (20 
seconds), 

2.  The  nitrogen  ejector  flow  was  turned  on  from  2  to  18  seconds, 

3 .  An  audio  signal  was  sent  to  the  videorecorder  at  0.5  seconds  as  a  time 
reference, 

4.  When  the  motor  chamber  pressure  reached  100  psia  a  1  second  delay  occurred, 
then  the  deflector  was  triggered  to  slide  out  of  the  way  (for  0.7  seconds)  and 
the  PDF  A  was  triggered  to  take  data  0.25  seconds  later,  and 

5.  At  the  termination  of  the  run  the  clock,  chamber  pressure,  three  diodes,  and 
deflector  position  data  were  output  to  a  separate  file. 


LABTECH  NOTEBOOK  Build-Time 


File  Run 


Figure  9.  LABTECH  Notebook  Data  Acquisition  and  Experiment  Conti'ol  Program 
for  PDPA  and  Three-Wavelength  Extinction  Measurements. 
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Keithley  MetraByte  hardware  provided  the  interfaces  between  the  computer  and 
the  equipment.  The  24-Bit  Parallel  Digital  VO  Board  (PIO-24)  and  the  Analog  &  Digital 
HO  Board  (DAS-16F)  were  internal  to  the  computer.  Externally,  the  PIO-24  was 
connected  to  a  24  Channel  Relay  Output  Board  (ERB-24).  The  DAS-16F  connected  to 
the  Screw  Terminal  Accessory  Board  (ST A- 16). 

G.  MALVERN  2600 

Ensemble  particle-sizing  was  done  using  the  Malvern  2600.  This  device  consists 
of  a  laser  transmitter  and  receiver  and  computer  for  control  and  data  acquisition.  A  9  mm 
diameter  He-Ne  laser  beam  with  wavelength  X  =  0.633  |xm  exits  the  transmitter  and  is 
scattered  (Fraunhofer  diffracted)  by  the  particles  based  on  their  size.  The  scattered  light 
then  passes  through  a  Fourier  transform  len  and  is  focused  on  32  semicircular  detector 
rings  to  measure  the  intensity  of  light  at  different  solid  angles  up  to  a  scattering  angle  0  of 
14°.  Due  to  the  physical  limitations  of  positioning  the  Malvern,  it  was  not  possible  to 
simultaneously  conduct  Malvern,  PDPA,  and  three-wavelength  extinction  measurements. 
Therefore,  the  Malvern  was  used  to  take  a  sample  through  the  plume  at  the  same 
longitudinal  position  6"  downstream  of  the  nozzle  as  the  main  PDPA  measurements.  It 
was  triggered  using  a  separate  LABTECH  Notebook  program  to  take  50  sweeps  of  data  in 
0.4  seconds,  starting  1  second  after  the  chamber  pressure  reached  100  psia.  [Ref.  17] 

H.  VIDEO  CAMERA 

A  video  camera  was  positioned  above  the  setup  and  framed  on  the  plume  region. 
The  signal  was  sent  live  into  the  control  room  where  it  was  monitored  and  recorded  on  a 
Sony  S-VHS  VCR.  Text  annotation  and  a  time  reference  were  displayed  on  the  screen, 
however  the  time  signal  was  not  synched  with  the  LABTECH  Notebook  clock. 

Therefore,  an  approximate  time  hack  was  sent  from  the  data  acquisition  program  in  the 
form  of  an  audio  signal  that  would  spike  the  VU  meter  on  the  VCR.  The  display  was 
used  to  measure  the  plume  dimensions,  determine  probe  tip  position  relative  to  the  plume 
centerline,  and  safely  observe  the  progress  of  the  run.  Although  1"  separated  holes  in  the 
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lab  table  were  typically  visible,  parallax  prevented  their  direct  use  in  measuring  motor 
nozzle  to  probe  tip  position. 
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IV.  EXPERIMENTAL  PROCEDURES 


Preparation  for  each  firing  required  several  steps  as  discussed  below.  Since  the 
igniters  required  extra  time  to  fabricate,  they  were  normally  made  in  advance  during 
non-firing  days.  Each  propellant  grain  was  cut,  bonded  to  the  motor  case  and  allowed  to 
cure  overnight.  On  the  day  of  the  firing  the  probe  and  side  tubes  were  assembled  and 
coarsely  aligned  on  their  mount  using  pencil  marks  of  the  final  previous  position  as  a 
guide.  The  argon  laser,  with  its  water  cooler  and  fan,  and  the  PDPA  were  then  energized 
to  conduct  fine  alignment.  Cigar  smoke,  deodorant  spray,  or  humidifier  water  vapor  were 
used  in  the  probe  to  see  the  beam  and  test  the  PDPA.  The  white  light  source  was  also 
turned  on  and  its  beam  checked  using  a  white  card  at  various  spots.  The  crucial  check 
was  a  clean,  collimated  beam  exiting  the  top  vertical  tube  of  the  probe.  Next,  the  diode 
box  was  coarsely  aligned  and  the  extension  tubes  were  put  in  place.  The  diode  box  was 
opened  and  the  filters  and  beam  splitters  cleaned  and  secured.  The  incoming  beam 
exiting  the  pinhole  was  positioned  on  the  center  of  each  diode.  The  cover  was  replaced 
and  the  diode  outputs  verified.  The  propellant  surface  was  seraped  to  provide  a  fresh 
surface  to  aid  in  ignition.  The  motor  was  assembled  without  the  igniter  and  mounted  at 
the  proper  longitudinal,  lateral,  and  vertical  position  in  front  of  the  probe  tip.  The  video 
camera  was  turned  on  and  adjusted.  The  probe  Nj  ejector  flow  valve  was  activated  to 
purge  the  line  and  then  connected  to  the  probe.  The  deflector  was  aligned  and  the 
activator  air  was  turned  on.  All  computers  were  booted  and  data  acquisition  software 
loaded.  The  videotape  recorder,  monitor,  and  screen  annotator  were  turned  on  and  set  up. 
The  white  light  beam  was  covered  temporarily  so  "no  light"  diode  readings  could  be 
taken.  Several  dry  runs  with  the  LABTECH  Notebook  program  were  performed  and  a 
manually  actuated  voltage  to  simulate  the  pressure  transducer  was  input  to  fully  check 
system  performance  and  data  acquisition.  Old  files  were  deleted  to  ensure  adequate  hard 
drive  space  for  the  0.5  MB  data  file.  Finally,  protective  shields  and  plastic  covers  secured 
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with  duct  tape  were  attached  to  protect  the  PDPA  from  exhaust  products.  The  igniter  was 
installed  and  battery/charger  hooked  up. 

The  firing  sequence  started  by  clearing  the  perimeter  area  of  personnel.  Next,  the 
warning  horn  was  sounded,  VCR  set  to  record,  PDPA  set  to  acquire  data  when  externally 
triggered,  and  the  LABTECH  Notebook  program  started.  The  final  event  was  a  short 
countdown  ''3-2-i-FIRE"  and  the  ignition  button  pressed  and  held  until  the  firing  was 
initiated  or  the  program  was  terminated  in  the  event  of  a  "no-fire". 

Initial  post-firing  procedures  consisted  of  stopping  PDPA  data  acquisition  and 
saving  the  data,  stopping  the  VCR,  and  securing  the  firing  circuit.  The  laser,  arc  lamp, 
diodes,  and  video  camera  were  turned  off.  The  battery  charger  was  unplugged  and  Nj/air 
valves  closed  for  the  ejector  and  deflector.  If  SEM  samples  were  to  be  taken,  then  they 
were  obtained  by  wiping  from  the  inside  of  the  back  of  the  probe  with  a  surgical  gloved 
finger  dipped  in  acetone.  Lenses,  mirrors,  and  other  glass  or  optics  were  wiped  with 
methanol  or  spray-cleaned  with  compressed  air.  All  exposed  areas  were  wiped  with 
acetone  to  try  to  prevent  corrosion  from  the  exhaust.  The  probe  and  all  tubes  as  well  as 
the  motor  were  disassembled  and  cleaned  with  water  or  acetone. 

A.  PHASE  1  FIRINGS 

Prior  to  the  operational  propellant  PS-1  firings,  several  test  firings  were  conducted 
with  AC- 14  to  check  the  setup  and  to  make  any  necessary  modifications.  Measurements 
using  the  PDPA  and  original  probe  were  made  and  drawings  made  of  the  required 
changes  and  additional  hardware.  Several  months  were  required  to  get  the  probe 
modified  (new  hole  cut  and  vertical  tubes  added)  and  the  PDPA  side  tubes  machined.  A 
collimation  tube  with  pinhole  and  lens  holders  was  manufactured  for  use  with  the  arc 
lamp.  Additional  probe  tips  (spares)  were  also  made  in  case  some  were  damaged  during 
firings. 

1.  AC-1410/24/94 

This  first  firing  was  to  check  whether  an  end-burning  grain  and  the  d,  =  0. 165" 
nozzle  would  yield  the  desired  250  psia  without  clogging.  It  also  was  an  opportunity  to 
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test  the  PDPA  through  the  newly  modified  probe  with  side  tubes.  The  probe  tip  was 
positioned  5.5"  (24  d^)  from  the  nozzle  and  1"  (4  d)  off  of  centerline.  Unfortunately,  the 
nozzle  clogged  with  molten  aluminum  and  silicon  forcing  the  chamber  pressure  up  to  600 
psia.  The  data  acquisition  program,  N2  ejector,  deflector,  and  PDPA  appeared  to  work  as 
desired,  although  the  trigger  times  needed  adjustment. 

2.  AC-1411/23/94 

The  second  firing  was  an  opportunity  to  try  a  different  grain  configuration  and 
nozzle  combination.  Two  grains,  one  end-  and  the  other  end-  and  radial-burning 
(described  in  Chapter  HI),  were  used  with  the  d,  =  0.239"  nozzle  to  ideally  achieve  350 
psia  for  about  1  second.  However,  a  somewhat  progressive  bum  resulted  in  a  chamber 
pressure  of  470  psia.  Clogging  was  still  a  problem.  The  probe  was  positioned  5.5"  (15 
d^  downstream  and  1.5"  (4  d^)  laterally.  The  PDPA  did  not  acquire  any  meaningful  data. 

3.  AC-1412/8/94 

This  was  the  first  run  with  the  three-wavelength  extinction  apparatus  in  addition 
to  the  PDPA.  A  small  diameter  tube  connected  the  top  vertical  tube  of  the  probe  to  the 
diode  box.  The  arc  lamp  was  positioned  on  a  stand  next  to  the  lab  table.  The  collimation 
tube  was  attached  to  the  condenser  housing  on  one  end  and  on  the  other  end  to  a  box 
containing  a  45°  mirror  beneath  the  lower  vertical  tube  of  the  probe.  It  was  discovered 
that  the  neutral  density  filters  and  beamsplitters  were  blocking  UV  light  (0.365  or  0.313 
)xm),  so  the  UV  filter  was  moved  to  the  lowest  diode  in  the  box  and  the  neutral  density 
filter  at  the  box  entrance  was  removed.  This  also  required  changing  neutral  density  filters 
in  front  of  the  other  diodes  until  suitable  output  levels  were  attained.  The  other  two-grain 
configuration  for  AC- 14  described  in  Chapter  HI  was  used  with  the  d^  =  0.268"  nozzle  to 
ideally  achieve  300  psia  for  1  second.  The  actual  pressure  peaked  at  420  psia  and  was 
regressive  burning  for  approximately  one  second.  Since  this  grain  and  nozzle 
combination  worked  quite  well  it  was  used  for  all  subsequent  runs  with  AC-14.  The 
diode  outputs  initially  started  to  drop  as  expected,  then  saturated,  then  decayed  to  a  new 
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full  light  value.  This  was  interpreted  as  a  vibration  or  shock  problem  resulting  from 
having  the  lamp,  collimator,  and  mirror  too  close  to  the  motor.  The  N2  ejector  line  had 
been  damaged  on  previous  runs  and  was  burned  off  this  time. 

4.  AC-141/18/95 

For  the  fourth  test  firing  several  changes  were  made.  A  hole  was  drilled  in  the 
table  and  the  white  light  source  was  moved  below  the  table.  An  adjustable  mirror  on  a 
post  with  a  magnetic  mount  was  added  below  the  table  and  replaced  the  mirror  in  the  box. 
It  was  also  determined  that  a  more  intense  and  clean  beam  resulted  without  the  secondary 
collimator,  so  it  was  no  longer  used.  The  light  was  collimated  instead  simply  by  the  lamp 
housing  lens  and  the  80"  distance  from  the  lamp  to  the  diodes.  The  probe  had  a  new  Nj 
ejector  connector  hole  drilled  on  the  opposite  side  for  better  protection.  The  pressure¬ 
time  trace  contained  a  1  second  neutral  bum  at  400  psia  and  the  diodes  responded 
properly  although  at  a  questionable  rate.  It  was  noted  that  the  diode  outputs  dropped 
during  Nj  ejector  activation,  probably  due  to  circulation  of  stray  particles.  The  vibration 
problem  appeared  to  be  solved.  A  hole  was  burned  in  the  side  funnel  which  connected 
the  probe  body  to  the  PDPA  detector. 

5.  AC-14  2/8/95 

The  side  fuimel  was  repaired  and  had  a  protective  plate  added  in  front.  The  diode 
response  question  was  handled  in  several  ways.  Unnecessary  amplifiers  that  had  a  1  Hz 
filter  were  disconnected.  The  diode  outputs  were  then  increased  by  adjusting  trim  pots 
and  swapping  out  neutral  density  filters.  The  sampling  rate  for  the  diodes  was  increased 
to  300  Hz  (this  also  applied  to  the  pressure  transducer  and  deflector  trigger).  Finally,  a 
light  chopper  was  used  to  produce  a  15  Hz  square  wave  and  the  diodes  had  no  difficulty 
measuring  the  alternating  light  conditions.  It  was  also  considered  that  a  sluggish  deflector 
could  have  caused  the  appearance  of  slow  diode  response.  Therefore,  the  video  replay 
was  used  to  verify  that  the  deflector  was  operating  as  commanded  as  accurately  as  could 
be  measured.  This  was  the  final  AC- 14  firing  because  everything  seemed  to  work 
properly.  The  acquisition  mn  time  was  decreased  from  30  to  20  seconds.  The  deflector 
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retraction  interval  was  increased  from  0.65  seconds  to  1  second.  The  probe  tip  was 
located  5.5"  (1 1  d)  from  the  nozzle  and  1"  (5  d^)  off  centerline.  The  N2  ejector  flow  was 
regulated  at  300  psia,  but  may  have  been  too  high  (non-isokinetic  flow)  as  evidenced  by 
the  diode  outputs  dropping  prior  to  the  deflector  retracting.  In  addition,  good 
transmittance  data  was  obtained  although  the  UV  diode  nearly  reached  T  —0.  This  was 
probably  due  to  too  large  a  regulated  Nj  pressure.  The  pressure-time  trace  showed  an 
average  of  350  psia  for  about  1  second.  Good  PDPA  data  was  obtained  (97  valid 

particles  with  a  d^2  - 

B.  PHASE  2  FIRINGS 

Since  the  AC- 14  firings  had  enabled  fixing  most  of  the  experimental  problems, 
the  operational  firings  with  PS-1  were  able  to  proceed  at  a  much  faster  pace.  The  new 
propellant  had  good  burning  characteristics,  ignited  easily,  and  did  not  produce  the  erratic 
"puffing"  of  exhaust  that  was  typically  observed  with  other  propellants.  It  also  was  easy 
to  cut  and  scrape.  However,  the  high  aluminum  content  of  the  propellant  produced  a  very 
difficult  to  remove  AI2O3  layer  in  the  nozzle. 

1.  PS-12/23/95 

As  with  AC- 14  an  end-burning  grain  was  initially  tried  for  simplicity  and 
assurance  of  neutral  burning.  The  d,  =  0.207"  nozzle  was  used  and  ideally  should  have 
resulted  in  a  chamber  pressure  of  224  psia  for  4  seconds.  However,  nozzle  clogging  was 
a  major  problem  and  a  highly  progressive  bum  resulted  with  a  peak  chamber  pressure  of 
430  psia.  The  probe,  PDPA,  and  three-wavelength  apparatus  were  not  used  since  this 
first  mn  was  for  plume  visualization  and  measurement.  The  clogging  problem  caused  the 
plume  to  be  irregular,  so  this  ran  was  not  as  useful  as  desired. 

2.  PS-12/25/95 

A  new  calibration  curve  was  generated  for  the  pressure  transducer  prior  to  this 
ran.  A  single  perforated  grain  and  the  larger  d,  =  0.239"  nozzle  were  selected  and  used 
for  the  rest  of  the  work.  Ideally  an  approximately  neutral  bum  near  300  psia  for  1.5 
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seconds  would  occur.  The  deflector  and  PDPA  triggering  were  changed  to  the  intervals 
0.25  -  1.25  and  0.45  -  1.25  seconds  after  reaching  100  psia  respectively.  The  N2  ejector 
regulated  pressure  was  decreased  to  100  psia.  The  motor  was  positioned  for  a 
nozzle-to-tip  distance  of  6"  (16  d^)  and  a  lateral  offset  of  1.2"  (3  The  diode  box 
pinhole  was  opened  slightly  creating  a  beam  with  a  diameter  of  about  0.29".  The  narrow 
top  tube  from  the  probe  to  the  diode  box  was  replaced  with  a  much  larger  1 "  O.D.  tube  to 
cut  down  on  internal  reflection,  improving  the  beam.  The  actual  chamber  pressure 
peaked  at  400  psia.  The  PDPA  did  not  measure  any  particles.  The  transmittances  looked 
good,  however  they  dropped  to  a  very  low  value.  Also,  because  the  transmittances  started 
dropping  early  it  was  felt  that  the  Nj  pressure  was  again  too  high.  From  the  video  replay 
it  appeared  that  the  tip  was  located  in  the  smoke  periphery  of  the  plume  as  desired.  The 
plume  width  was  estimated  to  be  3"  with  the  highly  visible  portion  about  1.5"  wide. 

3.  PS-12/27/95 

This  run  was  to  be  on  centerline.  Since  only  large  particles  {d  >  0.5  jxm)  were 
anticipated,  only  the  PDPA  was  used.  Concern  for  being  too  close  and  melting  the  tip 
drove  the  nozzle-to-tip  distance  up  to  12"  (33  d^).  To  reduce  the  tip's  exposure  time  to  the 
hot  plume  the  deflector  trigger  interval  was  shortened  to  0.25  -  0.95  seconds  after  100 
psia  chamber  pressure  was  reached.  The  actual  run  had  an  essentially  neutral  bum  at  a 
chamber  pressure  of  about  300  psia  for  2  seconds.  Unfortunately,  the  PDPA  only 
measured  8  valid  particles  with  a  d^^  =  12.0.  It  was  decided  that  the  three-wavelength 
measurements  should  always  be  taken  to  aid  in  evaluating  what  the  PDPA  measured.  The 
nozzle  throat  diameter  was  reduced  by  the  Al^Oj ,  and  was  measured  after  the  run  to  be 
0.15".  This  was  typical  of  all  the  subsequent  runs. 

4.  PS-13/1/95 

From  the  two  previous  mns  it  appeared  that  the  deflector  and  PDPA  were  being 
triggered  early,  so  the  intervals  were  again  changed  to  1  - 1.7  and  1.25  - 1.7  seconds  after 
100  psia  respectively.  The  nozzle-to-tip  distance  was  significantly  reduced  to  4"  (1 1 
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and  this  run  was  on  centerline.  No  PDPA  data  was  collected,  possibly  due  to  a  computer 
problem  after  an  automatic  laser  shutdown  occurred  after  an  overheat  condition  with  the 
water  cooler.  The  probe  tip  melted,  probably  0.3  seconds  after  deflector  retraction. 
Particles  were  collected  for  SEM  analysis. 

5.  PS-13/2/95 

The  nozzle-to-tip  distance  was  set  to  6"  (16  and  0.2"  (0.5  d)  laterally.  The 
chamber  pressure  peaked  near  380  psia  had  an  essentially  neutral  bum.  Transmittance 
data  from  this  and  the  previous  run  showed  that  the  ejector  flow  was  still  too  high.  Good 
PDPA  data  was  acquired  with  44  valid  particles.  Particles  were  collected  again  for  SEM 
analysis.  The  tip  melted,  probably  0.1  seconds  after  deflector  retraction. 

6.  PS-13/3/95 

For  this  run  the  nozzle-to-tip  distance  was  set  to  6"  (16  d)  and  0.4"  (0.5  d^ 
laterally.  The  ejector  flow  pressure  was  reduced  to  35  psia.  An  essentially  neutral  bum 
for  1.5  seconds  was  achieved  with  a  peak  pressure  of  360  psia.  The  PDPA  data  appeared 
to  be  correct.  However,  it  had  to  be  copied  by  hand  and  then  reconstmcted  due  to  a 
computer  hard  drive  problem.  Only  four  valid  particles  were  measured  with  a  d^2  -  2-3 
pm.  Transmittance  data  still  suggested  that  the  ejector  flow  rate  was  too  high. 

7.  PS-13/4/95 

This  mn  had  a  nozzle  to  tip  distance  of  6"  (16  d^  and  was  offset  0.6"  (1.6  d^ 
laterally.  The  ejector  flow  pressure  was  reduced  to  25  psia.  Due  to  rainy  weather  the 
motor  took  a  while  to  light-off,  but  did  achieve  a  neutral  bum  for  two  seconds  near  375 
psia.  The  PDPA  measured  only  one  valid  particle  at  6.4  pm,  and  the  transmittance  data 
looked  good. 

8.  PS-13/5/95 

In  order  to  conduct  a  background  radiation  check  and  try  an  open-air,  three- 
wavelength  measurement,  the  probe  with  side  and  vertical  tubes  as  well  as  the  PDPA 
were  not  used.  The  motor  was  set  up  for  a  6"  plume  length  and  a  0.75"  offset  from  the 
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light  beam  emanating  from  the  lab  table.  The  neutral  density  filters  and  diode  box 
pinhole  had  to  be  adjusted  to  prevent  the  diodes  from  saturating  with  the  increased 
radiance.  During  the  video  replay  a  faint  glow  of  the  beam  was  visible  where  it  came  out 
of  the  table  prior  to  ignition.  Once  the  motor  fired  the  beam  was  very  visible  as  it  was 
scattered  by  the  particles  in  the  plume.  The  plume  appeared  to  be  approximately  3" 
across  at  the  location  of  the  beam.  The  surprising  result  was  plume  emission  at  the  UV 
wavelength  X  =  0.313  )xm,  and  potentially  at  the  other  wavelengths. 

9.  PS-13/6/95 

Due  to  the  results  of  the  previous  run  a  light  chopper  was  added  to  the  setup 
between  the  light  source  and  mirror  and  run  at  a  speed  to  produce  an  approximately 
40  Hz  square  wave  of  alternating  light  and  darkness.  This  frequency  was  selected 
because  with  the  sampling  rate  of  300  Hz  it  was  the  highest  frequency  that  would  provide 
two  to  three  sample  points  near  the  top  and  bottom,  corresponding  to  light  and  no  light,  of 
the  wave.  The  motor  was  positioned  6"  (16  d)  from  the  probe  tip  on  centerline.  The 
bum  was  fairly  regressive  with  a  peak  of  500  psia.  The  PDPA  had  no  valid  particles  of 
192  attempts,  however  the  tip  melted  at  some  point  during  the  run.  Transmittance  data 
indicated  that,  inside  the  probe  plume,  emission  was  not  a  problem  since  throughout  the 
run  accurate  "no  light"  samples  were  recorded  at  the  bottom  of  the  square  wave. 

10.  PS-1  3/7/95 

The  final  firing  used  only  the  Malvern  particle  sizer.  The  Malvern  was  set  up  in  a 
separate  test  cell  with  a  100  mm  Fourier  transform  lens,  providing  a  vignetting  distance  of 
approximately  5"  in  front  of  the  receiver  lens.  A  laser  bandpass  filter  was  used  to 
eliminate  background  radiation.  The  motor  was  mounted  for  a  nozzle-to-beam  distance 
of  6"  with  the  plume  centerline  2.5"  from  the  receiver  lens.  A  separate  LABTECH 
Notebook  program  was  used  to  trigger  the  Malvern  1  second  after  100  psia  chamber 
pressure.  The  Malvern  was  initially  aligned  and  adjusted  for  as  low  a  background  as 
possible.  Then  it  was  set  in  a  program  mode,  awaiting  external  trigger  to  take  50  sweeps 
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in  0.4  seconds,  followed  by  a  10  second  delay,  then  another  50  sweep  sample  in  0.4 
seconds.  This  latter  sample  was  to  check  the  condition  of  the  windows  under  post-firing 
conditions.  The  video  camera  was  also  set  up  looking  down  on  the  plume.  The  motor 
did  not  have  much  of  a  neutral  bum  and  reached  a  peak  of  about  475  psia.  The  Malvern 
measurement  was  successful  and  occurred  near  the  peak  pressure.  Post-firing  inspection 
of  the  receiver  window  showed  that  it  was  partially  obscured  by  some  exhaust  particle 
residue.  The  video  replay  provided  an  excellent  opportunity  to  measure  the  plume 
dimensions  that  were  lost  on  the  first  motor  firing  with  PS-1.  The  highly  visible  portion 
of  the  plume  appeared  to  be  12"  long  and  1.5"  wide  with  the  widest  point  at  the  6" 
longitudinal  position.  The  outer  smoky  portion  of  the  plume  grew  from  about  the  3"  aft 
position  to  a  maximum  diameter  of  about  4"  near  the  6"  aft  position,  and  maintained  this 
width  to  the  edge  of  the  field  of  view. 
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V.  RESULTS  &  DISCUSSION 


This  chapter  examines  the  results  obtained  with  the  Malvern,  PDPA,  and  the 
three-wavelength  apparatus.  The  discussion  refers  to  the  PS-1  firings  on  2/25/95  and 
later.  Calculations  are  presented  to  correlate  the  data  by  different  means.  Table  4 
presents  the  PS-1  firings  that  were  used  for  data  reduction  (i.e.,  the  first  firing  has  not 
been  included).  The  d,  =  0.239"  nozzle  which  was  used  resulted  in  an  average  chamber 
pressure  fP, )  between  290  and  480  psia.  This  yielded  underexpanded  nozzle  flow  (exit 
pressure  between  30  and  40  psia  as  verified  with  the  MICROPEP  equilibrium,  adiabatic 
combustion/performance  code  run  in  the  off-design  mode  for  the  appropriate  expansion 
ratio). 


Firing 

Positiont 

PDPA 

3-X,  Ext. 

Malvern 

SEM 

Ave.  P, 
(psia) 

Comments 

IIIM^ 

6"/1.2" 

Yes* 

Yes 

No 

No 

380 

2/27/95 

12"/0.2" 

Yes 

No 

No 

No 

290 

3/1/95 

4'70" 

No 

Yes 

No 

Yes 

375 

Melted  tip 

6'70.2" 

Yes 

Yes 

No 

Yes 

360 

Melted  tip 

3/3/95 

6'70.4" 

Yes 

Yes 

No 

No 

325 

6’70.6" 

Yes 

Yes 

No 

No 

370 

6'70.75" 

No 

Yes 

No 

No 

400 

Bare  light 

6'70" 

Yes* 

Yes 

No 

No 

480 

Light  chopper, 
melted  tip 

3/7/95 

6"/NA 

No 

No 

Yes 

No 

470 

_  _  _ j  ^ 

Table  4.  Summary  of  PS-1  Firings  Used  for  Data  Reduction.  (*  -  An  attempt  was  made 
but  there  were  no  valid  particles,  t  ■  Axial  location/radial  displacement) 


A.  MOTOR 

As  discussed  in  the  previous  chapter  the  propellant  always  burned  in  a  smooth 
fashion.  Periodic  puffing  of  the  exhaust  from  shedding  of  Aip,  slag  was  not  observed. 
Formerly  used  propellants  have  shown  a  propensity  for  this  shedding  phenomenon  and 
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with  such  a  high  aluminum  content  in  PS-1  it  was  anticipated.  Chamber  pressures  varied 
between  250  and  500  psia  compared  to  the  near  300  psia  calculated  for  ideal  conditions, 
due  primarily  to  a  reduction  in  throat  area  from  AljOj  deposits. 

B.  EXPECTED  RESULTS 

Previous  experimentation  had  demonstrated  the  individual  usefulness  and 
limitations  of  the  PDPA,  the  Malvern,  and  the  multiple-wavelength  extinction  techniques 
for  the  solid  rocket  motor  environment.  The  Malvern  takes  an  across-the-plume 
measurement  of  "large"  particles  with  d  >  0.5  |im  and  provides  number  and  volume 
distributions  assuming  there  are  no  particles  smaller  than  0.5  |im.  The  PDPA  counts 
individual  particles  with  d  >  0.5  |xm  passing  through  its  probe  volume  and  provides 
number  and  velocity  distributions  as  well  as  representative  mean  diameters.  It  has  been 
demonstrated  to  accurately  measure  the  distribution  of  particles  with  d  >  0.5  |xm 
suspended  in  liquid.  Multiple-wavelength  extinction  measurements  have  also  been  used 
before  at  the  laboratory.  A  Mie  code  and  an  assumed  distribution  shape  (e.g.,  log-normal) 
permitted  the  particle  index  of  refraction  and  distribution  parameters  for  "small"  particles 
in  the  plume  edges  to  be  found.  Since  past  work  had  found  the  larger  particles  only 
toward  the  center  of  the  plume  and  the  smaller  particles  throughout  the  plume,  the  goal  of 
this  investigation  was  to  use  the  PDPA  to  obtain  "large"  particle  distributions  at  various 
radial  positions  across  the  plume  and  to  use  the  extinction  technique  at  the  plume  edge  to 
get  distribution  information  for  the  "small"  particles  which  then  could  be  applied  across 
the  whole  plume.  The  Malvern  was  to  be  used  to  provide  a  check  on  the  combined  radial 
PDPA  measurements.  In  order  to  accomplish  this  objective  experiments  were  conducted 
using  the  probe  at  radial  positions  of  0",  0.2",  0.4",  0.6",  0.75",  and  1.2"  from  the  plume 
centerline  and  a  final  measurement  was  made  across  the  plume. 

To  investigate  the  potential  for  taking  the  obtained  log-normal  distribution,  found 
for  the  small  particles  using  the  extinction  technique,  and  curve-fitting  it  to  a  PDPA 
distribution,  a  spreadsheet  was  created  in  LOTUS®  1-2-3®  Release  4  for  Windows™. 
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The  PDPA  distribution  was  manually  entered  for  all  the  possible  bin  diameters  and  bins 
for  additional  diameters  smaller  than  0.5  |xm  were  inserted  for  use  with  the  log-normal 
distribution.  The  equation  for  the  cumulative  log-normal  distribution,  Equation  (2.26), 
was  used  to  determine  the  number  of  particles  in  each  individual  bin  to  create  a  number 
distribution  in  histogram  form  for  comparison  to  that  of  the  PDPA.  An  overlapping 
portion  of  the  two  distributions  for  diameters  greater  than  0.5  jxm  was  used  for  the 
curve-fit.  A  Sum  of  the  Squared  Deviations  (SSD)  approach  was  used  to  obtain  a  single 
number  to  evaluate  the  curve-fit.  The  SOLVER  routine  in  1-2-3  was  used  to  adjust 
Gj,  and  //to  minimize  the  SSD.  Unfortunately,  even  with  constraints  the  routine  could 
not  accomplish  this.  It  could,  however,  work  with  a  single  parameter  if  the  other  two 
were  fixed.  By  manually  adjusting  the  parameters  the  problem  was  readily  apparent. 
Many  log-normal  distributions  can  provide  an  equally  good  fit.  Thus,  curve-fitting  will 
not  result  in  a  unique  solution. 

The  Mie  code  FORTRAN  program  at  NPS,  obtained  from  Cashdollar  [Ref.  12], 
was  formerly  only  available  on  the  main-frame  computer.  Therefore,  to  speed  up  the  data 
reduction  and  analysis,  it  was  copied  and  modified  for  use  on  a  personal  computer  with 
Microsoft  MS-Fortran  4.1.  To  simulate  experimental  conditions  in  the  1-2-3  spreadsheet 
and  examine  expected  transmittance  values  for  different  distributions  and  at  different 
wavelengths  "miescat3.exe"  (see  Appendix  C)  was  used  and  a  macro  was  written  for 

1-2-3  to  be  able  to  properly  import  the  output  file  (see  Appendix  D).  What  was  more 
important  was  also  to  be  able  to  examine  the  InT  or  Q  ratios  as  the  mass  of  "small" 
particles  was  varied  with  respect  to  the  mass  of  the  "large"  particles.  Finally,  for  data 

reduction  (i.e.,  finding  the  index  of  refraction  m,  and  g^  that  best  fits  the 
experimentally  derived  transmittances)  the  program  "miescat6.exe"  was  developed  (see 
Appendix  E). 

The  most  optimistic  result  desired  was  to  be  able  to  take  the  end-product, 
experimentally-derived  distributions  for  various  radial  positions  across  the  plume  and 
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conduct  a  mass  calculation  for  comparison  to  the  calculated  mass  that  should  be  present 
from  burning  the  propellant. 

C.  MALVERN 

The  Malvern  data  were  typical  of  past  observations  with  propellants  of  various 
composition.  Beam  steering  occurred  as  evidenced  by  the  increased  intensities  on  the 
first  six  detector  rings.  The  voltages  on  these  rings  were  deleted  before  the  size 
distribution  was  calculated.  Elimination  if  these  inner  six  rings  does  not  affect  the 
distribution  when  all  the  particles  have  diameters  less  than  the  expected  value  of 
approximately  30  pm.  In  agreement  with  previous  data  obtained  with  the  subscale  motor 
[Ref.  18],  the  Malvern  measured  ~25%  of  the  volume  in  particles  with  d  <  2.0  pm,  the 
remaining  -75%  of  the  volume  in  particles  with  d  =  A-5  pm,  and  a  d^  -  6  pm.  An 
obscuration  of  86%  (therefore  r=  0.14)  was  recorded,  with  beam-steering  present,  and  a 
^32  =  2.59  pm  was  determined. 

In  order  to  conduct  a  check  on  the  above  value  of  T,  Dobbins'  formula  (Equation 
2.30)  was  used.  First,  a  stagnation  enthalpy  energy  balance  equation  was  used  to  estimate 
the  ratio  of  the  mass  flow  rates  of  air  to  propellant  gas  in  the  plume  at  the  measurement 
location  for  use  in  calculating  the  mass  concentration  of  particles. 

^  P  gas  gas'll  gas  ^  p  air^ t  air  Q^Pgas  (5.1) 

where 

ihpg^  =  mass  flow  rate  of  propellant  gas 

Cpgas  =  specific  heat  for  propellant  gas  =  143 1  J/kg-K  from  MICROPEP 

Ttgas  =  stagnation  temperature  of  the  propellant  gas  =  3584  K  from 
MICROPEP 

ihair  =  mass  flow  rate  of  entrained  air 

^Pair  =  specific  heat  for  entrained  atmospheric  air  =  1004  J/kg-K 
=  stagnation  temperature  of  entrained  air  =  520  K 

Cpn^  specific  heat  of  the  mixture  of  gases  =1217  J/kg-K  from 
average  of  gas  and  air 

“  stagnation  temperature  of  the  mixture  of  gases  -1500  K  [Ref  19 
at  16 
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Dividing  by  the  mass  flow  rate  of  propellant  gas,  the  ratio  of  the  mass  flow  rates  was 
calculated.  The  mass  flow  rate  of  entrained  air  was  approximately  2.5  times  that  of  the 
propellant  gas.  Therefore,  the  mass  flow  rate  of  the  mixture  of  gases  was  3.5  times  the 
mass  flow  rate  of  the  propellant  gas  alone.  The  equation  used  for  C„  was 

^  ^PAIiOi  ^PAl20iPmix  _  fnpAl20iPmi^_ 

Cm  —  •  —  •  '  T?  T  ^ 

V^mix  ^mix^mix^  mix 

mix 


where 

Vmix=  volume  flow  rate  of  the  mixture  of  gases 
density  of  the  mixture  of  gases 
Pmij<  -  pressure  of  the  mixture  of  gases  ~  14.7  psia 
^mix-  constant  for  the  mixture  of  gases  ~  54  ft-lbf/lbm-°R 
temperature  of  the  mixture  of  gases  ~  1500  K 

First,  the  mass  flow  rate  of  the  propellant  had  to  be  calculated  using 

tkp  =  PpAbap"^  (5.3) 


where 

jhp  =  mass  flow  rate  of  the  propellant 

=  propellant  density  =  0.06665  Ibm/in^  from  MICROPEP 
A*  =  propellant  burning  surface  area  =  4.8  in*  average  during  steady  bum 
a  =  propellant  burning  rate  pressure  coefficient 
=  chamber  pressure  =  400  psia 
n  =  propellant  burning  rate  pressure  exponent 

Next,  the  mass  flow  rate  of  AlPj  was  found  by  using  the  mass  percentage  of  aluminum 
in  the  propellant  (20%)  times  the  ratio  (1.89)  of  the  molecular  weights  of  AljOj  (101.96 
g/mol)  to  aluminum  (2  x  26.98  g/mol)  times  the  mass  flow  rate  of  the  propellant  found 
using  Equation  (5.3).  Finally  all  values  were  plugged  into  Equation  (5.2)  to  find  C„=  1.4 
X  10'®  Ibm/in*.  Using  Equation  (2.30)  with  L  =  3",  Q  =  2,  p  =  Pai203=  0-138  Ibm/in*  at 
1500  K,  and  d^2  -  ^-59  |xm  a  transmittance  of  T  =  0.40  was  calculated.  This  is  much 

higher  than  the  Malvern  measured  value  of  0.14,  most  likely  due  to  beam  steering  as 
discussed  above.  Beam  steering  is  the  refraction  of  the  beam  due  to  the  thermal  gradient 
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in  the  path  of  the  beam.  This  causes  the  central  detector,  which  measures  the  obscuration, 
to  receive  less  intensity  than  it  should  and  the  innermost  detector  rings  to  receive  too 
large  an  intensity. 

Another  calculation  was  performed  to  compare  the  mass  in  the  observed  modes 
for  the  calculated  of  2.59  |im.  For  the  Malvern  volume  distribution  two  modes  were 
apparent,  near  1  and  4.5  |xm.  Therefore,  it  was  assumed  that  all  the  particles  were  either  1 
or  4.5  )i,m  and  the  number  ratio  of  1  to  4.5  p,m  particles  was  found  to  be  24.3  using: 

Nx  ^32(4.5)^ -(4.5)3 

A^4.5  (1)3-J32(1)2 


where  is  the  number  of  1  |j.m  particles  and  is  the  number  of  4.5  |im  particles. 
Then,  the  mass  percent  of  the  1  |im  particles  to  the  total  mass  of  the  particles  was 
calculated  as  21.1%  using 


Vi 

F total 


Ih. 

N4.5 


(1)^ 


Ni 

N4.5 


(1)3 +  (4.5)3 


(5.5) 


where  F,  is  the  volume  of  1  |im  particles  and  is  the  volume  of  all  the  particles.  Thus, 
the  Malvern  data  implied  that  most  of  the  mass  was  contained  in  particles  with  diameters 
greater  than  1  |xm.  This  is  evident  in  Figure  10  of  Appendix  G. 

D.  PDPA 

Six  attempts  were  made  to  obtain  measurements  at  various  radial  positions  in  the 
plume  with  the  PDPA,  from  the  centerline  out  to  1.2"  (3  d^).  Both  centerline  firings 
resulted  in  no  valid  particles.  On  both  of  these  runs  as  well  as  the  0.2"  offset  run  the 
probe  tip  was  melted.  Therefore,  flow  deflection  and/or  a  strong  detached  normal  shock 
could  have  caused  the  lack  of  "large"  particles  in  the  probe  by  breaking  them  up  into 
"small"  particles. 
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The  Malvern  number  and  volume  distributions  were  compared  to  an  equivalent 
across-the-plume  distribution  formed  by  scaling  the  0.2",  0.4",  and  0.6"  offset  PDPA 
distributions.  This  scaling  was  done  by  assuming  the  plume  was  axisymmetric  and  could 
be  divided  up  into  concentric  zones.  Within  each  zone  the  distribution  was  assumed  to  be 
constant.  Since  there  were  three  PDPA  distributions  that  could  be  used  three  zones  were 
formed  with  the  constraint  that  each  zone  could  contain  only  one  PDPA  distribution. 

This  forced  the  inner  zone  to  be  constrained  to  an  outside  radius  of  0.2  -  0.4".  Similarly, 
the  middle  zone  extended  from  the  outside  of  the  inner  zone  to  0.4  -  0.6"  and  the  outer 
zone  extended  from  the  outside  of  the  middle  zone  to  0.6  -  0.75".  Each  distribution  was 
then  multiplied  by  twice  the  difference  between  its  zone's  inside  and  outside  radii.  This 
was  the  scale  factor.  Next,  the  scaled  distributions  were  combined  and  the  total  number 
of  particles  found.  This  total  was  then  divided  into  each  number  value  and  multiplied  by 
100  to  obtain  number  %  at  a  given  diameter.  Finally,  the  number  %'s  were  added 
together  according  to  the  bin  diameters  used  by  the  Malvern  and  plotted  with  the  Malvern 
distribution.  An  additional  step  was  required  for  the  volume  distribution.  Each  scaled 
number  was  multiplied  by  (P.  Then  when  combined,  the  total  volume  of  particles  was 
found  rather  than  number.  The  PDPA  and  Malvern  plots  of  the  number  and  volume 
distributions  (see  Figures  10  - 13  in  Appendix  G)  agreed  that  there  was  a  mode  at  <  1  p,m 
and  another  mode  between  4  and  5  fim.  The  PDPA  volume  distribution  highlights 
individual  larger  particles  that  it  measured,  whereas  individual  particles  are  essentially 
lost  by  the  Malvern  as  it  looks  at  an  ensemble  of  particles. 

The  PDPA  distributions  also  showed  that  there  were  almost  no  "large"  particles  (d 
>  0.5  |xm)  outside  of  0.6"  off  the  centerline.  This  was  in  agreement  with  the  visual 
appearance  of  the  plume.  The  brighter  plume  core  area  appeared  to  have  approximately  a 
1"  diameter. 

There  are  several  ways  to  improve  the  future  performance  of  the  PDPA  when  used 
in  the  plume  enviroiunent.  First  and  foremost,  if  the  probe  tip  were  cooled  or  made  of  a 
more  thermally  resistant  material  the  tip  erosion  problem  could  be  minimized.  A  shorter 


49 


time  for  tip  exposure  would  also  help  in  this  regard.  The  ejector  flow  needs  to  be 
optimized  to  prevent  sucking  in  too  many  of  the  small  particles  since  this  biases  the  data 
toward  smaller  particles.  It  would  be  very  difficult  to  obtain  isokinetic  sampling  because 
of  the  required  short  explosive  time.  The  probe  needs  to  be  redesigned  with  a  smaller 
frontal  area,  probably  circular,  and  with  a  diameter  no  larger  than  the  present  maximum 
tip  diameter.  This  would  cut  down  on  the  plume  obstruction  and  remove  some  of  the 
excess  internal  volume  that  presently  causes  recirculation  within  the  probe.  If  fiber  optics 
were  used,  the  large  side  tubes  could  also  be  removed,  further  improving  the  plume 
obstruction  problem.  With  subscale  motors  it  may  also  be  possible  to  not  use  the  probe, 
allowing  the  PDPA  to  traverse  the  plume  cross-section. 

E.  THREE-WAVELENGTH  EXTINCTION 

Transmittances  were  calculated  for  seven  firings  and  used  to  obtain  "best  fit" 
values  of  index  of  refraction  m,  a^,  and  for  a  log-normal  distribution.  The  following 
formula  was  used: 


Vp-Vo 

Vf-Vo 


(5.6) 


where  is  the  diode  output  in  volts  for  particle  obscured  light,  Vp  is  the  diode  output  for 
no  light,  and  V^is  the  diode  output  for  unobstructed  full  light.  The  no-light  reading  was 
made  prior  to  firing  by  covering  the  hole  in  the  bottom  of  the  lab  table.  The  full  light 
value  was  found  by  taking  an  average  of  the  first  two  seconds  of  the  run  (i.e.,  before 
ejector  gas  activation  which  typically  reduced  the  diode  reading  shghtly).  The  particle 
obstructed  reading  was  taken  over  a  flat  minimum  during  steady-state  burning. 

Once  the  three  transmittances  were  calculated  for  a  firing  they  were  inserted  into 
the  data  input  file  "miescat6.dat"  and  iteration  intervals  were  chosen  for  m  and  a^.  The 
data  reduction  file,  "miescat6.exe",  was  run  and  provided  the  best  for  o^oh 

combination  of  m  and  as  well  as  a  calculation  of  the  degree-of-fit.  This  value  was  the 
Sum  of  the  Squared  Deviations  (SSD)  between  InT  ratios  and  Q  ratios.  The  closer  SSD 
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approached  zero  the  better  the  fit.  Table  5  presents  all  of  this  information.  The  SSD 
values  are  presented  as  integers  rather  than  the  eight  digits  to  the  right  of  the  decimal  that 
they  really  were  (i.e.,  237  is  really  237E-08).  The  3/5/95  firing  was  different  from  the 
others  because  the  probe  was  not  used  and  the  plume  was  definitely  emitting  at  0.313  |xm. 
Therefore,  a  transmittance  value  was  not  calculated  at  this  wavelength  and  a  fit  could  not 
be  obtained. 

It  was  discovered  that  the  process  used  to  iterate  toward  the  "best-fit"  values  was 
convergent  for  and  but  not  for  m.  Therefore,  since  a  reasonable  value  for  m  of 
AljOj  in  the  plume  is  1.70  -  Oi,  this  value  was  used  as  a  starting  point  to  find  a  reasonable 


Firing 

Position 

^313 

T 

^546 

T 

m 

^32 

SSD 

2/25/95 

6" / 1.2" 

0.0557 

0.0533 

0.0958 

1.70-0i 

1.82 

0.6 

237 

3/1/95 

4"  /  0" 

0.2891 

0.3304 

0.4905 

1.70-0i 

1.7 

0.44 

1,708 

3/2/95 

6"  /  0.2" 

0.4957 

0.5054 

0.5666 

1.70-0i 

2.1 

0.69 

925 

3/3/95 

6"  /  0.4" 

0.3026 

0.425 

0.5567 

1.70-0i 

3.24 

0.36 

292 

3/4/95 

6"/ 0.6" 

0.2518 

0.3315 

0.4316 

1.70-0i 

3.38 

0.58 

253 

3/5/95 

6"/ 0.75" 

N/A 

0.794 

0.7184 

N/A 

N/A 

N/A 

N/A 

3/6/95 

6"  /  0" 

0.1711 

0.2268 

0.3007 

1.70-0i 

3.9 

0.84 

207 

Table  5.  Three- Wavelength  Extinction  Results. 


value  for  and  d^2-  Then,  if  a  better  fit  was  desired,  m  had  to  be  varied  over  its  entire 
possible  range  in  0.01  increments  as  was  stepped  in  its  smallest  increment  of  0.01 
about  the  reasonable  value  and  the  minimum  SSD  set  chosen.  This  entire  somewhat 
lengthy  process  was  not  done,  however  some  excellent  fits  were  found.  For  the  2/25/95 
firing  m  =  1.79  -  Oi,  =  1.84,  =  0.55  resulted  in  an  SSD  =  lE-08.  For  the  3/4/95 

firing  m  =  1.85  -  Oi,  =  3.38,  t/jj  =  0.50  resulted  in  an  SSD  =  2E-08.  The  large  values  of 

Og  and  d22  relative  to  the  0.5  |J.m  particles  is  really  indicative  of  a  significant  amount  of 
the  mass  in  "large"  particles.  It  should  be  noted  that  only  two  independent  ratios  were 
available  using  three  wavelengths.  Kim  et  al.  [Ref.  20]  used  six  wavelengths  for  five 
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independent  ratios  and  determined  a  ~  1.50  and  ~  0150  p.m  in  the  edge  of  the 
plume.  Thus,  if  further  work  is  to  be  done,  it  should  use  as  many  wavelengths  as  possible 
with  a  large  difference  between  X  ^  and  k 

An  uncertainty  analysis  was  conducted  on  the  transmittances,  then  applied  to  a 
single  transmittance  and  propagated  through  the  fit  process  to  determine  how  much  the 
distribution  parameters  changed.  Using  Holman's  [Ref.  21]  approach  the  transmittance 
uncertainty  was  found  with 


Wt- 


y-pvr\ 
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Vp-yf 

[Vf-Vof 


r^'vo 


Vo-Vp 

[Vf-VoY 


^2 


:WVf 


(5.7) 


where  w^,p,  w^^,  and  Wjy.are  the  uncertainties  in  the  transmittance,  particle  obscured 
diode  output,  no-light  diode  output,  and  full-light  diode  output,  respectively.  Using  the 
values  from  the  3/4/95  data  set  w^,  and  were  set  equal  to  0.0397,  the  1  a 
deviation  on  the  546  diode  between  0  and  2  seconds  elapsed  time.  This  resulted  in  a 
transmittance  uncertainty  of  0.008974.  This  was  added  to  to  create  a  perturbed  set, 
then  the  fit  process  was  conducted  using  m  =  1.70  -  Oi.  This  produced  a  =  5.85, 
d-^2  =  0.71,  and  an  SSD  =  8E-08.  If  this  is  compared  to  the  original  fit  in  Table  5  it  is 
obvious  that  these  values  correspond  to  two  unique  distributions.  Therefore,  the  multiple 
wavelength  extinction  technique  is  highly  sensitive  to  any  experimental  uncertainties. 

The  3/5/95  bare-light  data  (no  probe)  were  used  to  compare  the  mass  of  small 
particles  to  large  particles.  Since  the  light  beam  was  0.75"  from  the  centerline  of  the 
plume  and  appeared  in  the  video  replay  to  be  in  the  smoky  portion  of  the  exhaust,  it  was 
assumed  to  be  in  the  small-particle  region.  Equation  (2.30)  was  used  to  determine  C„. 
The  plume  was  estimated  to  be  3"  in  diameter,  so  L  was  calculated  to  be  2.60".  The 

distribution  parameters  found  using  the  fit  routine  with  the  2/25/95  data  taken  in  the 
plume  edge  were  used  with  the  Mie  code  to  find  Q.  However,  the  transmittance  values 
were  not  used  because  the  high  probe  purge  rate  resulted  in  a  high  concentration  of  small 

particles  in  the  probe  (note  low  values  of  T).  Using  the  3/5/95  value,  a  C„  of  1.2E-07 
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Ibm/in^  was  calculated.  When  compared  to  the  total  C„  of  1.4E-06,  calculated  above, 
this  showed  that  the  small  particles  accounted  for  only  8%  of  the  total  mass.  Similarly, 
using  the  3/5/95  value  a  C„  of  6.4E-08  Ibm/in^  was  calculated  which  corresponds  to  a 
small  particle  mass  of  only  4%  of  the  total.  When  the  mass  distribution  is  plotted  for  the 
small  particles  (see  Figure  18  in  Appendix  G)  notice  that  only  30%  of  the  "small  particle" 
mass  is  in  particles  with  diameter  less  than  0.5  p-m.  Therefore,  the  mass  %  of  particles 
with  d  <  0.5  pm  is  less  than  3%.  Thus,  accurate  measurement  of  the  size  distribution  in 
the  plume  should  be  possible  using  only  the  PDF  A.  This  greatly  simplifies  the 
diagnostics  required  to  make  the  needed  plume  particle  size  measurements. 

F.  SEM  ANALYSIS 

Several  samples  of  the  particles  coating  the  inside  rear  of  the  probe  were  taken  for 
analysis  under  a  scanning  electron  microscope  (SEM).  First,  the  particles  were 
sedimented  in  a  beaker  with  acetone  to  try  to  remove  any  other  matter.  Next,  they  were 
scraped  from  the  bottom  and  sprinkled  on  carbon  painted  mounts.  Finally,  the  mounts 
were  placed  in  a  vacuum  container  where  carbon  fibers  were  vaporized  to  carbon  coat  the 
samples.  This  was  to  help  with  the  sample  charging  problem.  The  samples  were 
examined  qualitatively  under  the  SEM  for  an  idea  of  the  general  size  distribution.  See  the 
photographs  in  Figure  19  of  Appendix  G.  It  was  noted  that  there  were  many  1  to  2  and  4 
to  6  p,m  particles  and  few  smaller  or  larger  particles.  This  is  consistent  with  the  measured 
distribution,  using  the  Malvern  and  the  PDPA ,  as  well  as  the  various  mass  calculations. 
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VI.  CONCLUSIONS 


The  PDPA  and  Malvern  distributions  agreed  in  the  observed  modes  near  1  and 
4.5  |im.  SEM  pictures  of  collected  particles  were  in  good  qualitative  agreement  with  the 
Malvern  Sauter  mean  diameter  =  2.59  p.m.  Data  analysis  indicates  that  less  than  3% 
of  the  total  mass  of  the  particles  was  contained  in  particles  with  diameter  d  <  0.5  |xm. 

This  means  that  the  PDPA,  which  can  typically  measure  particles  down  to  a  minimum 
diameter  of  0.5  p,m  with  a  dynamic  range  (d^:d^„)  of  50:1,  can  be  used  by  itself  to 
determine  the  particle  size  distribution.  Although  multiple-wavelength  measurements 
were  not  required,  they  were  found  to  be  very  sensitive  to  inaccuracies  in  the  measured 
transmittances.  The  modified  particle  collection  probe  allowed  simultaneous  PDPA  and 
multiple-wavelength  measurements.  However,  the  design  should  be  improved  in  several 
ways.  The  frontal  area  of  the  probe  should  be  reduced  to  lessen  plume  obstruction.  One 
option  is  to  use  a  small  cylindrical  cross-section  and  use  fiber  optics  with  the  PDPA, 
thereby  removing  the  large  side  tubes  which  also  obstruct  the  plume.  In  addition,  the  tip 
exposure  time  should  be  reduced  or  a  more  thermally  resistant  material  chosen  for  the  tip. 
Further  studies  should  be  conducted  to  check  shock  formation  in  and  around  the  tip.  A 
"rake"  of  PDPA  miniature  fiber  optic  probes  may  be  a  reasonable  design  to  investigate  for 
taking  simultaneous  radial  measurements.  Traversing  such  a  rake  across  the  plume  is 
einother  option. 
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APPENDIX  A.  MICROPEP  OUTPUT 


The  following  pages  contain  a  printout  from  the  equilibrium,  adiabatic 
combustion/performance  code  MICROPEP.  This  particular  run  was  for  propellant  PS-1 
with  a  chamber  pressure  of  400  psia  and  a  nozzle  area  ratio  of  2.33  operating  off-design 
for  an  ambient  sea  level  pressure  of  14.7  psia.  The  area  ratio  pertains  to  the  d,  =  0.239" 
nozzle  used  for  most  of  the  experiments  with  PS-1.  Note  that  the  exit  pressure  is  about 
40  psia,  indicating  underexpanded  flow. 
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Feb.  1994  -  Modified  by  A.  McAtee 

Naval  Postgraduate  School,  Monterey,  CA 

****  NEWPEP  **** 

*  PS-1  Ae/At=2.33  * 


INGREDIENT 

MASS 

HF 

DENSITY 

COMPOSITION 

(gm) 

(kcal/kg) 

(kg/in**3 ) 

AMMONIUM  PERCHLORATE  (AP) 

70.00 

-602.0 

1948.6650 

ICL 

4H 

IN 

ALUMINUM  (PURE  CRYSTALINE) 

20.00 

.0 

2701.5590 

lAL 

HTPB  (SINCLAIR) 

10.00 

13.0 

918.9728 

103H 

73C 

10 

VOLUME  PERCENT  OF  INGREDIENTS  (IN  ORDER) 
66.2684  13.6572  20.0744 


THE  PROPELLANT  DENSITY  IS  .06665  LB/CU-IN  OR  1.8448  GM/CC 
THE  EQUIVALENCE  RATIO  IS  1.3260 

NUMBER  OF  GRAM  ATOMS  OF  EACH  ELEMENT  PRESENT  IN  INGREDIENTS 

3.416525  H  .732471  C  .595760  N  2.393073  0 

.741290  AL  .595760  CL 


****************************CHAMBER  RESULTS  FOLLOW***************************** 

TEMP  PRESSURE  ENTHALPY  ENTROPY  CP/CV  SGAMMA  Pi/ni 

(K)  (MPa/ATM/PSI)  (kJ/kg)  (kJ/kg-K)  (MPa/kmol) 

3584.7  2.758/  27.22/  400.00  -1757.6980  9.398  1.1620  1.1212  84385.270 

DAMPED  AND  UNDAMPED  SPEED  OF  SOUND=  848.508  AND  1063.794  m/sec 

SPECIFIC  HEAT  (MOLAR)  OF  GAS  AND  TOTAL=39670 . 840  54036.950  J/kmol-K 

NUMBER  MOLS  GAS  AND  CONDENSED=  3.2682  .3392 


*=liguid 

.84621 

,  &=solid) 
H2 

.68459 

CO 

.48272 

H20 

.45869 

HCl 

.33919 

A1203* 

.29499 

N2 

.22731 

H 

.08121 

Cl 

.06428 

HO 

.04771 

C02 

.03125 

AlCl 

.01142 

AlOCl 

9.90E-03 

0 

5.82E-03 

A1C12 

5.60E-03 

NO 

4.92E-03 

A1H02 

3.22E-03 

AlO 

2.99E-03 

AlHO 

2.21E-03 

02 

1.46E-03 

A1 

5.45E-04 

A120 

3.80E-04 

A1C13 

1.56E-04 

AlH 

1.12E-04 

C12 

8.25E-05 

OCl 

7.47E-05 

A1202 

6.17E-05 

N 

5.99E-05 

CHO 

4.18E-05 

A102 

4.07E-05 

HOCl 

2.93E-05 

COCl 

1.82E-05 

NH 

1.81E-05 

NH3 

1.81E-05 

HO  2 

1.75E-05 

NH2 

8.93E-06 

CNH 

7.17E-06 

NHO 

4.92E-06 

AlHO 

2.68E-06 

CH20 

l.OOE-06 

CNO 

THE  MOLECULAR  WEIGHT  OF  THE  MIXTURE  IS  27.721  gm/mole 
THE  GAS  CONSTANT  IS  299.93  JAg-K 

TOTAL  HEAT  CONTENT  (298  REF)  =  6284.328  kJ/kg 

SENSIBLE  HEAT  CONTENT  (298  REF)=  5721.090  kJAU 

************************* ***exHAUST  results  FOLLOW***************************** 

TEMP  pressure  ENTHALPY  ENTROPY  CP/CV  SGAMMA  Pi/ni 

(K)  (MPa/ATM/PSI)  (kJ/kg)  (kJ/kg-K)  (MPa/kmol) 

2917.2  .276/  2.72/  40.02  -3746.0840  9.398  1.1556  1.1201  8780.586 

DAMPED  AND  UNDAMPED  SPEED  OF  SOUND=  732.222  AND  938.454  m/sec 
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SPECIFIC  HEAT  (MOLAR)  OF  GAS  AND  TOTAL=39437 . 610  55329.410  J/kmol-K 

NUMBER  MOLS  GAS  AND  CONDENSED=  3.1425  .3642 


(*=liquid,&=solid) 


.88420 

H2 

.67999 

CO 

.53227 

HCl 

.48500 

H20 

.36417 

A1203* 

.29725 

N2 

.12299 

H 

.05240 

C02 

5.05E-02 

Cl 

2.15E-02 

HO 

6.85E-03 

AlCl 

3.03E-03 

AlOCl 

1.82E-03 

0 

1.27E-03 

A1C12 

1.20E-03 

NO 

8.34E-04 

A1H02 

4.60E-04 

AlHO 

4.01E-04 

02 

2.01E-04 

AlO 

1.33E~04 

A1C13 

9.42E-05 

A1 

3.24E-05 

Cl  2 

1.53E-05 

A120 

8.73E-06 

OCl 

6.01E-06 

CHO 

5.81E-06 

AlH 

5.79E-06 

HOCl 

4.70E-06 

N 

3.41E-06 

COCl 

2.96E-06 

NH3 

1.86E-06 

A1202 

1.42E-06 

NH2 

1.15E-06 

A102 

1.09E-06 

HO  2 

1.03E-06 

NH 

l.OOE-06 

CNO 

THE  MOLECULAR  WEIGHT  OF  THE  MIXTURE 

IS  28.517 

gm/mole 

THE  GAS  CONSTANT  IS  291.56  J/kg-K 

TOTAL  HEAT  CONTENT  (298  REF)  =  4788.339  kJ/kg 

SENSIBLE  HEAT  CONTENT  (298  REF)=  4456.853  kJ/kg 

******************************************************************************* 


An  exact  method  for  determining  throat  conditions  was  used 

The  frozen  &  shifting  STATE  gammas  for  the  throat  are:  1.1600  1.1199 

GAMMA  NU  shown  below  is  the  gamma  for  the  chamber  to  throat  PROCESS. 


**********PERF0RMANCE:  FROZEN  ON  FIRST  LINE,  SHIFTING  ON  SECOND  LINE,********* 

**************SHIFTING  TO  THROAT/FROZEN  AFTERWARDS  ON  THIRD  LINE************** 


SPECIFIC 

GAMMA 

T* 

P* 

C* 

ISP* 

Pe 

D-ISP 

A*/m 

Te 

Cf 

IMPULSE 

NU 

(gm-s/ 

(cm* *2/ 

(sec) 

(K) 

(MPa) 

(m/s) 

(sec) 

(MPa) 

cm**3) 

kg/s) 

(K) 

222.8 

1.1630 

3314. 

1.577 

1539.2 

.276 

10192.710 

5.5805 

2590. 

1.2718 

226.8 

1.1202 

3416. 

1.601 

1560.4 

195.8 

.297 

10238.540 

5.6571 

2937. 

1.2602 

224.9 

1.1202 

3416. 

1.601 

1560.4 

195.8 

.297 

10140.370 

5.6571 

2702. 

1.2482 

FROZEN  & 

SHIFTING  KINETIC 

ENERGY  1 

OF  EXHAUST  674708.  781129.  m**2/s**2 
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ia-141  so  sNtm 

13.143  lOO  SHKETS 
33-144  300SHCETS 


?DPA  ecvi^ 
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'd3n  09  7 
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TVP  —  pLOrV\^  PRoB^ 


APPENDIX  C.  MIE  CODE  "MIESCAT3.FOR  " 


On  the  following  pages  is  the  Fortran  program  to  calculate  extinction  coefficients 
Q,  average  extinction  coefficients  Q,  and  ratios  of  average  extinction  coefficients  for  the 
inputs  given  in  "iniescat3.dat".  The  original  program  was  run  on  the  main-frame 

computer  at  NPS.  Since  the  Combustion  Research  Lab  (CRL)  did  not  have  a  network 
connection  to  the  main-frame  the  program  was  copied  and  modified  to  run  on  an  IBM 
compatible  PC  using  the  Microsoft®  MS®  Fortran  4. 1  compiler  at  the  CRL. 

The  data  block  in  the  original  program  followed  the  end  subroutine.  Since  the 
source  code  had  to  be  compiled  and  linked  prior  to  execution,  it  was  inconvenient  to  have 
to  repeat  these  steps  each  time  the  inputs  were  changed.  Therefore,  the  data  block  was 
removed  and  put  in  the  separate  file  "miescat3.dat"  and  the  read  statements  modified  in 
the  source  code. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


BY  KENNETH  CASHDOLLAR,  1976 
REVISED  BY  THOMAS  WELDON,  JULY,  1977 
BASED  ON  PROGRAM  MICOEF  BY  C  D  LITTON  AND  Z  J  FINK 
REVISED  BY  LT  JOHN  MANSER,  USN  1994-1995  (UPDATED  FOR 
MS-FORTRAN  4 . 1 , DATA  BLOCK  WITH  WAVES  PUT  INTO 
EXTERNAL  FILE) 


"MIESCAT3.EXE"  CALCULATES  EXTINCTION  COEFFICIENTS  AND  RATIOS  FOR 
LIGHT  TRANSMISSION  THROUGH  DUST  PARTICLES  WITH  COMPLEX  REFRACTIVE 
INDEX  AND  A  LOG-NORMAL  SIZE  DISTRIBUTION 


**********  "MIESCAT3.DAT”  INPUT  DATA  FILE  FORINT  ****************** 


VARIABLE  NAMES  AS  IN  PROGRAI^ 

ROW  1 . WAV(l) ,WAV(2) ,WAV(3) 

ROW  2 . NDATA 

ROW  3 . RDUST,RMED 

ROW  4 . L2,L3,DS,DELDS 

ROW  5 . ND, SIGMA, SD32,SDL32 

FORMATS  AND  VARIABLE  DESCRIPTION 

ROW  1,...F6.4  (WAVELENGTH  1),2X,F6.4  (WAVELENGTH  2),2X, 

F6.4  (WAVELENGTH  3) 

ROW  2.... 13  (DATA  SET  IDENTIFIER,  IF  EQUAL  TO  999  SUBSEQUENT 
DATA  IS  NOT  PROCESSED) . 

ROW  3....2F6.3  (COMPLEX  REFRACTIVE  INDEX  OF  DUST  PARTICLES), 
2X,F6.3  (REFRACTIVE  INDEX  OF  THE  MEDIUM). 

ROW  4.... I 3  (L2,  THE  NUMBER  OF  INCREMENTS  TAKEN  AT  A  SET  DELTA 
)(THE  NliMBER  OF  TIMES  DELTA  IS  DOUBLED,  L3),F6.3 
(STAR  DIAMETER) ,F6. 3 (INITIAL  VALUE  OF  DELTA). 

ROW  5 - 13 (NUMBER  OF  INCREMENTS  TAKEN  FOR  LOG-NORM) , F4 . 1 ( SIGM 

F5 . 2 ( STARTING  DIAMETER ) , F5 . 2 ( DELTA ) . 


c*********** 

c 

c 

c 

c 


COMPLEX  R , ROUST , DI , DELDI , RX , RINDX , DIS , DELIS 
COMPLEX  COEFR,DPR(1000) ,XJR(1000) ,DJR(1000) ,PR(1000) 
DIMENSION  XY(IOOO) ,WAV(3) 

DIMENSION  DP(IOOO) ,XK(1000) ,XJ(1000) ,DJ(1000) ,P(1000) 
DIMENSION  RATl(lOO) ,RAT2(100) ,RAT3(100) 

DIMENSION  PLOTDy(lOO) 

COMPLEX  A(IOOO) ,B(1000) ,C(1000) 

COMPLEX  ALPH,XH(1000) ,DH(1000) ,E(1000) ,DE(1000) 
COMMON/ADQ/DQ (500,2) 

COMMON  ALLQ(I00,4,3) 

COMMON/KPLT/RINDX ( 3 ) , KAVL ( 3 ) 
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oonoo  o  on  oooo  on  o  noon  oo  oooo 


COMMON  DIAM32(100) 


IPEN  =  4 

CALL  SHERPA( 'FILENAME' , 'M' , IPEN) 

OPEN  (UNIT=1,  FILE='MIESCAT3 .DAT' ) 
READ  (1,103)WAV(1) ,WAV(2) ,WAV(3) 
103  FORMAT(F6.4,2X,F6.4,2X,F6.4) 

DO  8  NDATA=1,10 
READ  (1,105)NDAT 
105  F0RMAT(I3) 

IF(NDAT.EQ.999)  GO  TO  9 


RDUST  IS  COMPLEX  REFRACTIVE  INDEX  OF  DUST  PARTICLES,  WITH  MINUS 
RMED  IS  REAL  REFRACTIVE  INDEX  OF  MEDIUM 
READ  (1,32)RDUST,RMED 
32  FORMAT(2F6.3,2X,F6.3) 

RX  IS  RELATIVE  REFRACTIVE  INDEX,  (N-IK)/M 

RX=RDUST/RMED 

IMAGINARY  PART  OF  R  CHANGED  TO  POSITIVE  TO  CONFORM  TO  EXPANSION 

FORMULA  IN  PROGRAM 

R=CONJG(RX) 

RINDX(1)=RX 

RINDX(2)=RX 

RINDX(3)=RX 

L2*L3  =  NUMBER  OF  SIZE  INCREMENTS 

INITIAL  DIAMETER  IS  DS,  INITIAL  VALUE  OF  DELTA  IS  DELDS 

READ  ( 1, 34 )L2,L3,DS, DELDS 
34  FORMAT(2I3,2F6,3) 

INITIAL  VALUE  OF  COMPLEX  DIAMETER  IS  DIS 
DIS=CMPLX(0.0,DS) 

DELIS=CMPLX (0.0, DELDS ) 

READ  (1,12)ND, SIGMA, SD32,SDL32 
12  FORMAT(I3,F5.2,2F5.2) 

WAVE  IS  WAVELENGTH  OF  LIGHT  IN  MICRONS 
WAV(1)=0.3650 
WAV(2)=0.5461 
WAV(3)=0.8520 
INITIAL  VALUES  FOR  ARRAY 
DO  51  K=l,3 
DO  52  J=l,2 
ALLQ(1 , J,K)=0.0 
ALLQ( 1 , J+2 ,K)=0 .0 
DO  53  1=  2,100 
ALLQ(I,J,K)=100.0 
ALLQ(I,J+2,K)=1.0 
53  CONTINUE 
52  CONTINUE 
51  CONTINUE 

DO  54  1=1,100 
DIAM32(I)=0.0 
RAT1(I)=0.0 
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RAT2(I)=0.0 
RAT3(I)=0.0 
PLOTDY(I)=0.0 
54  CONTINUE 
DO  4  NI=1,3 
C 

PRINT  106,NDATA 

106  FORMAT(/'  DATA  SET  NUMBER', 13) 

PRINT  35 

35  FORMAT('  MIE  COEFFICIENTS  FOR  MONODISPERSE  DUST  PARTICLES') 

C 

PRINT  41,RDUST,RMED,RX 

41  FORMATC'  RDUST=' ,2F7.4,/'  RMED=' , F7 . 4 , /'  REL,  REF.  INDEX= ' , 2F7 . 4 ) 
PRINT  104, R 

104  FORMAT('  R==',2F7.4) 

WAVE=WAV(NI) 

PRINT  30, WAVE 

30  FORMAT('  WAVELENGTH  =',  F7 . 4 , 2X , 'MICRONS '/ ) 

PRINT  25 

25  FORMAT(6X, 'A' ,12X, ' D' , 13X , 'QABS ' , IIX , 'QEXT'/) 

WAVL(NI)=WAVE*10000 . 0 
C 

c 

C  INITIAL  VALUES  SET 

C 

C  MAX  VALUE  FOR  D  IS  ( 1000*WAVE) / ( 1 . 3*PI*RMED )  =  244*WAVE/RMED 

C  D  IS  REAL  DIAMETER  IN  MICRONS, DELD  IS  INCREMENT 

C  DI=I*D  IS  COMPLEX  DIAMETER,  DELDI  IS  INCREMENT 

C  X  IS  REAL  SIZE  PARAMETER  X=  PI*DIAMETER*RMED/WAVELENGTH 

C  ALPH  IS  COMPLEX  SIZE  PARAMETER,  ALPH-I*X 

D=DS 

DELD=DELDS 

DI=DIS 

DELDI=DELIS 


M=1 

DO  3  J=1,L3 

DO  2  1=1, L2 

X=3.14159*D/WAVE*RMED 

ALPH=3 . 14159*DI/WAVE*RMED 

XJ{1)=  SIN(X)/X 

XJ(2)=  SIN(X)/(X**2)-  COS(X)/X 

XJR(1)=CSIN(R*X)/(R*X) 

XJR(2)=CSIN(R*X)/( (R*X)**2)-CCOS(R*X)/(R*X) 

XH ( 1 ) =CEXP ( ALPH ) /ALPH 

XH(2)=( (-CEXP(ALPH) ) /X) -( ( 0 . , 1 . ) *CEXP( ALPH)/ ( X**2 ) ) 
DJ(1)=  COS(X)/X-  SIN(X)/(X**2) 
DJR(l)=CCOS(R*X)/(R*X)-CSIN(R*X)/( (R*X)**2) 
DH(l)=CEXP(ALPH)/X+( (0. , 1 . ) *CEXP( ALPH) /( X**2 ) ) 
QEXT=0 . 0 
QSCA=0 . 0 

LI  =  ITERATIONS  FOR  SPHERICAL  FUNCTIONS,  L1=1,3*X 
L1=1.3*X 

IF(Ll.LE.lO)  Ll=10 
DO  1  L=1,L1 
CO=(FLOAT(L) )/(2*L+l) 

COE= ( FLOAT ( L+1 ) ) / ( 2  *L+1 ) 
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C0EF=(FL0AT(2*L+1) )/X 
C0EFR={FL0AT(2*L+1) )/(R*X) 

XJ ( L+2 ) =COEF*XJ ( L+1 ) -XJ ( L) 

X  JR  ( L+2 )  =COEFR*X  JR  ( L+1 )  -X  JR  ( L ) 

XH  ( L+2 )  =COEF*XH  { L+1 )  -XH  ( L ) 

DJ(L+l)=C0*XJ(L)-C0E*XJ(L+2) 

DJR ( L+1 ) =CO*XJR( L) -COE*XJR( L+2 ) 

DH ( L+1 ) =CO*XH ( L ) -COE*XH ( L+2 ) 

P(L)=X*XJ(L+1) 

PR(L)=(R*X)*XJR(L+1) 

DP(L)=XJ(L+1)+X*DJ(L+1) 

DPR(L)=XJR(L+1)+(R*X)*DJR(L+1) 

E(L)=X*XH(L+1) 

A(L)=( (DPr|l)*P(L) )-(R*PR(L)*DP(L) ) )/( (DPR(L)*E(L) )-(R*PR(L)*DE(L 

^B(L)=( (R*DPR(L)*P(L) )-(PR(L)*DP(L) ) )/( (R*DPR(L)*E(L) )-(PR(L)*DE(L 
1))) 

C(L)=A(L)+B(L) 

GAMMA=CABS ( A ( L ) ) 

BETA=CABS(B(L)  ) 

XK(L)=(FLOAT(2*L+l))*(2./(X**2))*((GAMMA**2)+(BETA**2)) 

XY ( L)=( FLOAT ( 2*L+1) )*(2./(X**2) )*(REAL(C(L) ) ) 

QEXT=QEXT+XY(L) 

QSCA=QSCA+XK(L) 

1  CONTINUE 

QABS=QEXT-QSCA 

DQ(M,1)=D 

DQ(M,2)=QEXT 

PRINT  22,X,DQ(M,1) ,QABS,DQ(M,2) 

22  FORMAT(lX,F8.3,6X,F7.3,9X,F6.3,9X,F6.3) 

D  AND  QEXT  FOR  MONODISPERSE  PUT  INTO  PLOTTING  ARRAY 
EVERY  FOURTH  POINT  CALCULATED  IS  PUT  INTO  ARRAY-  ALLQ 


XM=M 

AREL=XM/4 . 0 

MINT=M/4 

XMINT=MINT 

IF  (AREL.NE.XMINT)  GO  TO  60 
IF  (D.GT.1.20)  GO  TO  60 

C  NUMBER  OF  POINTS  FOR  PLOTS  1  AND  3  IS  NPTSl 

NPTS1=MINT 
ALLQ(MINT+1,1,NI)=D 
ALLQ ( MINT+1 , 3 , NI ) =QEXT 
ALLQ ( MINT+1 , 1 , NI ) =D 
ALLQ (MINT+1 , 3 ,NI )=QEXT 
60  CONTINUE 

C  INCREMENT  SIZE  PARAMETERS  D=DIAMETER 

D=D+DELD 
DI=DI+DELDI 
M=M+1 

2  CONTINUE 

C  PRINT  36 

36  FORMAT  (IH  ) 

DELD=2.0*DELD 

DELDI=2.0*DELDI 

3  CONTINUE 
PRINT  106,NDATA 
PRINT  235 
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235  FORMAT('  MIE  COEFFICIENTS  FOR  LOG-NORMAL  DUST  DISTRIBUTION') 

PRINT  2 36, SIGMA 

236  FORMAT  ('  SIGMA=' , F4 . 2 ) 

PRINT  237,RX 

237  FORMAT('  RELATIVE  REFRACTIVE  INDEX=' , 2F7 . 4 ) 

PRINT  30, WAVE 

CALL  QAVG(M,NI,SIGMA,NPTS2,SD32,SDL32,ND) 

K1=ND+1 

DO  499  N=1,K1 

PLOTDY(N)=ALLQ(N,4,NI) 

499  CONTINUE 

IF(NI.EQ.2)  GO  TO  502 
IF(NI.EQ.3)  GO  TO  503 

501  CALL  PLOTD(DIAM32, PLOTDY, 41, .FALSE. , 'LINLIN' , 'NOLEGEND' , 

+  '  ', 'Diameter,  microns' , 'Extinction  Coefficient') 

GO  TO  4 

502  CALL  PLOTD (DIAM3 2, PLOTDY, 41, .FALSE. , 'LINLIN' , 'NOLEGEND' , 

+  '  ', 'Diameter,  microns' , 'Extinction  Coefficient') 

GO  TO  4 

503  CALL  PLOTD (DIAM3 2, PLOTDY, 41, .TRUE.  , 'LINLIN' , 'NOLEGEND' , 

4*  '  ', 'Diameter,  microns ',' Extinction  Coefficient') 

4  CONTINUE 

PRINT  106,NDATA 

7  CONTINUE 
PRINT  500 

500  FORMAT('  QBAR  RATIOS  FOR  3  WAVELENGTHS') 

PRINT  2 3 6, SIGMA 

PRINT  237,RX 

PRINT  510,WAV(3) , WAV( 1 ) , WAV( 3 ) ,WAV(2) ,WAV(2) ,WAV(1) 

510  FORMAT (/2X,'  D32  ' , 3 ( 5X , 'QBAR( ' , F6 . 4 , ' ) /QBAR( ' , F6 . 4 , ' ) ' ) / ) 

K=ND+1 

DO  600  N=2,K 

RAT1(N)=ALLQ(N,4,3)/ALLQ(N,4,1) 

RAT2(N)=ALLQ(N,4,3)/ALLQ(N,4,2) 

RAT3(N)=ALLQ(N,4,2)/ALLQ(N,4,1) 

PRINT  650,DIAM32(N) ,RAT1(N) ,RAT2(N) ,RAT3(N) 

600  CONTINUE 

650  FORMAT(lX,F6.3,13X,F7.4,23X,F7.4,23X,F7.4) 

PRINT  350,NDATA 

350  FORMAT  (/'**************************  END  OF  DATA  SET  NUMBER', 13/) 
CALL  PLOTD ( DIAM3  2(2), RATI ( 2 ) , 40 , . FALSE . , ' LINLIN ' , ' NOLEGEND ' , 

+  '  ', 'Diameter,  microns ', 'Extinction  Coefficient') 

CALL  PLOTD(DIAM32(2) ,RAT2(2) ,40, .FALSE. , 'LINLIN' , 'NOLEGEND' , 

+  '  ', 'Diameter,  microns' , 'Extinction  Coefficient') 

CALL  PLOTD(DIAM32(2) ,RAT3(2) ,40, .TRUE.  , 'LINLIN' , 'NOLEGEND' , 

+  '  ', 'Diameter,  microns' , 'Extinction  Coefficient') 

8  CONTINUE 

9  CONTINUE 

CALL  DONEPL 
END 


SUBROUTINE  QAVG ( M , NI , SIGMA , NPTS2 , SD3 2 , SDL3 2 , ND ) 

BY  KENNETH  CASHDOLLAR,  1976 

QAVG  CALCULATES  D32  AND  QBAR  FOR  A  GIVEN  QEXT  VS  D  AND 
A  GIVEN  LOG-NORMAL  SIZE  DISTRIBUTION 


INPUT  DATA  =  ND,SIGMA,D32,DEL32 
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c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 


COMMON  ALLQ(100,4,3) 

COMMON/ADQ/DQ (500,2) 

COMMON  DIAM32(100) 

DQ(I,1)=D,  DQ(I,2)=QEXT 

NI=WAVELENGTH  OR  REFRACTIVE  INDEX  CHANGE 
ND  IS  NUMBER  OF  D32,QBAR  CALCULATIONS 
SIGMA  IS  THE  GEOMETRIC  STANDARD  DEVIATION 
D32  IS  THE  SURFACE  WEIGHTED  MEAN  DIAMETER 
DELS 2  IS  INCREMENT  FOR  D32 

DIST  IS  SURFACE  WEIGHTED  LOG-NORMAL  SIZE  DISTRIBUTION 
M  IS  THE  MAXIMUM  NUMBER  OF  POINTS  IN  Q-ARRAY  FOR  MONODISPERSE 

PRINT  240 

240  FORMAT(3X, 'D32' ,8X, 'DMIN' ,7X, 'DMAX' ,5X, 'NMBR' ,6X, 'DlOX' ,7X, 'D32X' 
1,7X, 'QBAR'/) 

D32=SD32 

DEL32=SDL32 

NUMBER  OF  POINTS  FOR  PLOTS  2  AND  4  IS  NPTS2 

NPTS2=ND+1 

DO  17  NN=1,ND 

DMIN=D32/SIGMA**3. 

DMAX=D32*SIGMA**2.5 

AS=ALOG(SIGMA)**2. 

DG=EXP ( ALOG (D32)-2.5*AS) 

QSUM=0 . 0 
DSUM=0 . 0 
XSUM=0 . 0 
DHOLD=0 . 0 
DNSUM=0 . 0 
XNSUM=0 . 0 
NMBR=0 
MM=M-2 

DO  13  N=1,MM 

IF  (DQ(N,1) .LT.DMIN)  GO  TO  13 
IF  (DQ(N,1) .GT.DMAX)  GO  TO  15 
IF  (NMBR.GT.O)  GO  TO  242 
DMINZ=DQ(N,1) 

242  CONTINUE 

DD=DQ(N+1,1)-DQ(N,1) 

AA=(ALOG(DQ(N+l,l)/DG) )**2/(2.0*AS) 

DIST=EXP ( -AA) *DQ ( N+1 , 1 ) /SQRT (6.28  3  2*AS ) 

QSUM=QSUM+ ( DIST*DQ ( N+1 , 2 ) +DHOLD*DQ ( N , 2 ) ) *DD/2 . 0 
DSUM=DSUM+(DIST*DQ(N+1 ,l)+DHOLD*DQ(N, 1) )*DD/2 .0 

XSUM=XSUM+ ( DIST+DHOLD ) *DD/2 . 0 

XNSUM=XNSUM+ ( DIST/DQ ( N+1 , 1 ) *  *  2+DHOLD/DQ (N,1)**2)/2.0*DD 
DNSUM=DNSUM+ ( DlST/DQ ( N+1 , 1 ) +DHOLD/DQ (N,1))/2.0*DD 
XHOLD=DQ(N+l,l) 

DHOLD=DIST 
NMBR=NMBR+1 
13  CONTINUE 
15  CONTINUE 

QBAR=QSUM/XSUM 
D32X=DSUM/XSUM 
D10X=DNSUM/XNSUM 
ALLQ(NN+1,2,NI)=D32X 
ALLQ ( NN+1 , 4 , NI ) =QBAR 
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DIAM32(NN+1)=D32 

PRINT  241,D32,DMINZ,XHOLD,NMBR,D10X,D32X,QBAR 
241  F0RMAT(1X,3(F6.3,5X) , 13 , 5X , 2 ( F6 . 3 , 5X) , F7 . 4 , 5X) 
D32=D32+DEL32 
17  CONTINUE 
RETURN 
END 
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APPENDIX  D.  MACRO  FOR  LOTUS  1-2-3  RELEASE  4 


On  the  following  pages  is  a  macro  for  Lotus®  1-2-3®  Release  4  for  Windows'*'^. 
When  the  output  data  file  from  "miescat3.exe"  was  imported  into  a  1-2-3  spreadsheet  all 
the  data  would  enter  as  a  single  column  due  to  text  blocks  between  sets  of  data. 

Therefore,  data  parsing  was  necessary  to  split  the  data  and  enter  it  into  the  appropriate 
number  of  columns.  Since  this  operation  needed  to  be  performed  every  time  the 
"niiescat3.exe"  data  was  to  be  used  it  was  advantageous  to  create  a  macro  to  automate  this 
task.  The  macro  was  recorded  as  the  data  parsing  was  done  manually  for  one  entire  file, 
then  the  macro  was  edited  until  it  was  streamlined  and  all  the  bugs  were  fixed.  It  is 
included  here  specifically  for  those  who  must  also  work  with  the  "miescat3.exe"  output. 

It  has  been  subsequently  discovered  that  by  using  the  File  Open,  [miescat3.exe 
output  filename].  Combine,  Formatted  Text  combination  that  the  data  will  be  imported 
properly.  However,  all  of  the  text  blocks  including  data  labels  will  be  missing. 

Therefore,  the  use  of  this  macro  may  be  the  best  option. 


73 


PL 


5 


CD 


TJ 

c. 

at 


X  4 

Os 


t-  O 

5  ^ 

<  < 

^  o 

2>  Q 

<  < 

5<o< 

5<a< 


?  5 


ui  o  P  o 

WuJliiUJ 

_J  -J  _J 

,  LU  LU  LU 
Q.  (O  Q  W 


^  ^ 

^  ^ 

O  CM 
ISJ  5  > 

X<  CM 

Z  CM  S 
O  5  < 

<<2 

<  UJ 
CO  a  CO 


c_  * 

X  > 

O  r  ; 


S  fsl- 


r*  ■< 

*  o 

-  -  lO 
a>  CM 
:?;  > 


5;  ~  - 

<<  * 

?! 

§|3|: 
I?  5<‘¥- 
tiiSGtiil 

Hi  O  ID  LU  1 
—1  ^  -J  _l  . 
tli  <  LU  LU  t 
Q  Q.  CO  O  ( 


_i 

<  <  < 

^  <  ;-■ 

s*i 

N  CM  > 

c23<q2 

>  ^  I—  LU  ^ 
Ui  tU  O  K  LU 

^  WuJUJ  w 

fc?iiiiijS 

CO  Q.  CO  Q  CL 


Q 

< 

ili 

;s< 

:  t£  < 

l  UJ  LU 
1  —I  _J 
I  UJ  UJ 
)  Q  CO 


A 

A 

A 

$. 

« 

I 


if 

‘^2  =  .  CM 

L|i 

|55- 

z  —  CO 
Nl  lO  > 

QC  <  < 

o<  ^ 


11 

r*  < 

s  ® 

d>  S 

CO  > 
IC  — 

<  < 
<  - 


>  ^  I—  LU 
LU  LJJ  O  h“ 
_»  O  LU  LU 
>-  QC  ^  — 1 
H-  <  UJ  UJ 
CO  0.  CO  Q 


5; 

LU  O 

<  LU 
CL  CO 


-  A 

S-SSf; 

^  tr  'T-  0 

gigs 

<  DC  < 

^  •  o  < 

i<o^ 

LU  I—  T  ^ 
FO  LU 
LU  UJ  U  CO 
— I  ^  >>  C]C 
UJ  lU  h-  < 
Q  CO  CO  0. 


Si 


r' 

A 


<  ^  < 
S|8§: 
<oS<( 
<oc5<i 

J —  LU  ^  L 

UJ  LU  CO  Ui  i 
-J  -J  DC  ^  . 
LU  U  <  LU  L 


r  Qi  ^  < 

:  O  <  ^ 

7  Z  Q>  CO  §  < 

co2<oS 

-  T  ^  LU 

>LU  WOFL 

■  j  CO  LU  LU  c 


<  DC  ^  > 


ICO  a 


LU  UJ  « 
CO  cu  c 


iS§S?S 

<5  J<oc 

I  h-  >  ^  H“  LU 
I  Olu 

I  LU  LU  LU 

I  -J  >-  DC  -i  ^ 
I  LU  I-  <  UJ  UJ 
I  CO  CO  CL  CO  O 


^-cd 


S|2 

<K  < 

h-  LU  F 
O  F  O 
LU  LU  LU 
^  — 1  ^ 
UJ  UU  UJ 
CO  Q  CO 


O  ®  p 
N  P  > 
DC  <  < 
O  <  ^ 

5p 

^<< 
lU^O 
^  CO  UJ 

^iiu 

CO  CL  CO 


^^25’ 


<^< 

m 

«S3  5- 


A  r\  rs 

r»  «  r« 

r*  x^ 

^il-|  I  I  g'ill: 

SLi?  a3  <<< 

®  z  •  -  IC 
CM  ^  o  2: 

CO  p  F  go 

(3  oQ  > 


fcfcL. 
S§5? 
<  <  < 
lU  LU  LU 
>  >  > 

pis 

ssss 

LU  D  3  3 

F  O  O  O 

uJfff 

JQQQ 
!)  LU  LU  LU 


APPENDIX  E.  MIE  CODE  "MIESCAT6.FOR" 


To  simplify  the  data  reduction  process  (i.e.,  finding  the  best  combination  of 
complex  index  of  refraction  m,  geometric  standard  deviation  and  Sauter  mean 
diameter  the  measured  transmittances)  the  Mie  Code  in  Appendix  B  was  further 
modified.  Each  modification  resulted  in  the  program  getting  assigned  a  new  number  such 
as  "miescat4.for".  Since  the  input  data  file  and  the  compiled  executable  version  usually 
also  changed  they  were  both  given  the  new  number.  This  particular  program  allows 
iteration  starting  values,  step  sizes,  and  the  number  of  loops  to  be  input  in  the  data  file  for 
m  and  a^.  The  program  outputs  the  current  values  of  m  and  and  the  best  fit  as  well 

as  the  Sum  of  the  Squared  Deviations  (SSD)  value  for  the  fit  between  the  InT  ratios  and 
the  Q  ratios. 
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BY  KENNETH  CASHDOLLAR,  1976 
REVISED  BY  THOMAS  WELDON,  JULY,  1977 
BASED  ON  PROGRAM  MICOEF  BY  C  D  LITTON  AND  Z  J  FINK 
REVISED  BY  LT  JOHN  MANSER,  USN  1994-1995  (PLOT 
FUNCTIONS  REMOVED,  UPDATED  FOR  MS-FORTRAN  4.1, 
DATA  BLOCK  WITH  WAVES  &  T  PUT  INTO  EXTERNAL  FILE, 
DATA  REDUCTION  ITERATIONS  ADDED) 


"MIESCAT6.EXE"  CALCULATES  EXTINCTION  COEFFICIENTS  AND  RATIOS  FOR 
LIGHT  TRANSMISSION  THROUGH  DUST  PARTICLES  WITH  COMPLEX  REFRACTIVE 
INDEX  AND  A  LOG-NORMAL  SIZE  DISTRIBUTION  AND  OUTPUTS  THE  SUM  OF 
THE  SQUARED  DEVIATIONS  BETWEEN  THE  QBAR  AND  LN(T)  RATIOS 


****************  INPUT  DATA  FILE  **MIESCAT6 .  DAT*’  FORMAT  ************* 


VARIABLE  NAMES  AS  IN  PROGRAM 

ROW  1. . . .T1,T2,T3 

ROW  2. . . .WAV(l) ,WAV{2) ,WAV(3) 

ROW  3. . . .NUMR,STEPR 

ROW  4 . . . . NUMS , STEPS 

ROW  5. . . .RDUST,RMED 

ROW  6, . . .L2,L3,DS,DELDS 

ROW  7. . . .ND,SIGMA,SD32,SDL32 

FORMATS  AND  VARIABLE  DESCRIPTION 

ROW  ONE - F6. 4 (TRANSMITTANCE  FOR  WAVELENGTH  1 ) , 2X, F6 . 4 (TRANSM 

ITTANCE  FOR  WAVELENGTH  2 ) , 2X , F6 . 4 ( TRANSMITTANCE  FOR 
WAVELENGTH  3) 

ROW  TWO - F6. 4 (WAVELENGTH  1 ), 2X , F6 . 4 (WAVELENGTH  2),2X,F6.4 

(WAVELENGTH  3)  (ALL  IN  MICRONS) 

ROW  THREE. .13 (NUMBER  OF  ITERATIONS  ON  RDUST ) , 2X , F4 . 3 ( STEP  SIZE 
OF  ITERATIONS). 

ROW  FOUR. . .13 (NUMBER  OF  ITERATIONS  ON  SIGMA) , 2X , F4 . 3 ( STEP  SIZE 
OF  ITERATIONS). 

ROW  FIVE. . .2F6. 3 (COMPLEX  REFRACTIVE  INDEX  OF  DUST  PARTICLES), 
2X,F6. 3 (REFRACTIVE  INDEX  OF  THE  MEDIUM). 

ROW  SIX. . . .2I3(L2,THE  NUMBER  OF  INCREMENTS  TAKEN  AT  A  SET  DELTA) 
(L3,  THE  NUMBER  OF  TIMES  DELTA  IS  DOUBLED) , 2F6 . 3 
(STAR  DIAMETER  AND  INITIAL  VALUE  OF  DELTA). 

ROW  SEVEN. .13 (NUMBER  OF  INCREMENTS  TAKEN  FOR  LOG-NORM) , F4 . 2 
(SIGMA) ,2F5. 2 (STARTING  DIAMETER  AND  DELTA). 


C* ********************************** *********************************** 


COMPLEX  R , RDUST , RDUSTl , DI , DELDI , RX , RINDX , DIS , DELIS 
COMPLEX  COEFR,DPR(1000) ,XJR(1000) ,DJR(1000) ,PR(1000) 
DIMENSION  XY(IOOO) ,WAV(3) 

DIMENSION  DP(IOOO) ,XK(1000) , XJ( 1000 ), DJ( 1000 ) ,P(1000) 
DIMENSION  RATI (100) ,RAT2(100) ,RAT3(100) 

DIMENSION  RAT4 ( 100 ), RATS (100) ,RAT6(100) ,SSD(100) 
COMPLEX  A(IOOO) ,B(1000) ,C(1000) 
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COMPLEX  ALPH,XH(1000) ,DH(1000) , E( 1000 ) ,DE ( 1000 ) 
COMMON/ ADQ/DQ (500,2) 

COMMON  ALLQ(100,4,3) 

COMMON /KPLT/RINDX ( 3 ) , WAVL ( 3 ) 

COMMON  DIAM32(100) 


OPEN  (UNIT=1,  FILE=' MIESCAT6.DAT') 

READ  (1,103)T1,T2,T3 

READ  {1,103)WAV(1) ,WAV(2) ,WAV(3) 

103  FORMAT(F6.4,2X,F6.4,2X,F6.4) 

READ  (1,104)NUMR,STEPR 
READ  (1, 104 )NUMS, STEPS 

104  FORMAT(I3,2X,F4.3) 

DO  8  NDATA=1,10 
READ  {1,105)NDAT 
105  FORMAT(I3) 

IF(NDAT.EQ.999)  GO  TO  9 
NDAT=001 

RDUST  IS  COMPLEX  REFRACTIVE  INDEX  OF  DUST  PARTICLES,  WITH  MINUS 
RMED  IS  REAL  REFRACTIVE  INDEX  OF  MEDIUM 
READ  (1,32)RDUST,RMED 
32  FORMAT(2F6.3,2X,F6.3) 

READ  (1,34)L2,L3,DS,DELDS 
34  FORMAT (2 13, 2F6. 3) 

READ  (1, 12 )ND, SIGMA, SD32,SDL32 
12  FORMAT(I3,F5.2,2F5.2) 

RDUST1=RDUST 
PRINT  18 

18  FORMAT ('  EXP.  TRANSMITTANCES : ' ) 

PRINT  20,WAV(1) ,T1,WAV(2) ,T2,WAV(3) ,T3 

20  FORMAT(3('  T ( ' , F6 . 4 , ' )=' / F6 . 4 ) /) 

PRINT  21 

21  FORMAT ('  BEST  LOG-NORMAL  FIT  VALUES:') 

PRINT  22 

22  FORMAT(8X, 'M' ,11X, 'SIGMAG' ,7X, 'D32' ,10X, 'SSD' ) 

DO  9  INDEX1=1,NUMS 

DO  8  INDEX2=1,NUMR 

RX  IS  RELATIVE  REFRACTIVE  INDEX,  (N-IK)/M 

RX=RDUST/RMED 

IMAGINARY  PART  OF  R  CHANGED  TO  POSITIVE  TO  CONFORM  TO  EXPANSION 

FORMULA  IN  PROGRAM 

R=CONJG(RX) 

RINDX(1)=RX 

RINDX(2)=RX 

RINDX(3)=RX 

L2*L3  =  NUMBER  OF  SIZE  INCREMENTS 

INITIAL  DIAMETER  IS  DS,  INITIAL  VALUE  OF  DELTA  IS  DELDS 


INITIAL  VALUE  OF  COMPLEX  DIAMETER  IS  DIS 
DIS=CMPLX(0.0,DS) 

DELIS=CMPLX (0.0, DELDS ) 

INITIAL  VALUES  FOR  ARRAY 
DO  51  K=l,3 
DO  52  J=l,2 
ALLQ(1,J,K)=0.0 
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ALLQ(1, J+2,K)=0.0 
DO  53  1=  2,100 
ALLQ(I, J,K)=100,0 
ALLQ(I, J+2,K)=1.0 

53  CONTINUE 
52  CONTINUE 
51  CONTINUE 

DO  54  1=1,100 

DIAM32(I)=0.0 

RAT1(I)=0,0 

RAT2(I)=0.0 

RAT3(I)=0.0 

54  CONTINUE 
DO  4  NI=1,3 

PRINT  106,NDATA 

106  FORMAT(/'  DATA  SET  NUMBER', 13) 

PRINT  35 

35  FORMAT('  MIE  COEFFICIENTS  FOR  MONODISPERSE  DUST  PARTICLES') 

PRINT  41,RDUST,RMED,RX 

41  FORMAT('  RDUST=' ,2F7,4,/'  RMED= ' , F7 . 4 , /'  REL.  REF.  INDEX=' , 2F7 . 4 ) 
PRINT  104, R 

104  FORMAT('  R=',2F7.4) 

WAVE=WAV(NI) 

PRINT  30, WAVE 

30  FORMAT('  WAVELENGTH  =' ,  F7 . 4 , 2X, 'MICRONS '/) 

PRINT  25 

25  FORMAT (6X, 'A' ,12X, 'D' ,13X, 'QABS' ,11X, 'QEXT'/) 

WAVL ( NI ) =WAVE* 1 0000.0 


INITIAL  VALUES  SET 

MAX  VALUE  FOR  D  IS  ( 1000*WAVE)/( 1 . 3*PI*RMED)  =  244*WAVE/RMED 

D  IS  REAL  DIAMETER  IN  MI CRONS, DELD  IS  INCREMENT 

DI=I*D  IS  COMPLEX  DIAMETER,  DELDI  IS  INCREMENT 

X  IS  REAL  SIZE  PARAMETER  X=  PI *DIAMETER*RME D/WAVELENGTH 

ALPH  IS  COMPLEX  SIZE  PARAMETER,  ALPH=I*X 

D=DS 

DELD=DELDS 

DI=DIS 

DELDI=DELIS 


M=1 

DO  3  J=1,L3 

DO  2  1=1, L2 

X=3 . 14159*D/WAVE*RMED 

ALPH=3 . 14159*DI/WAVE*RMED 

XJ(1)=  SIN(X)/X 

XJ(2)=  SIN(X)/(X**2)-  COS(X)/X 
XJR(1)=CSIN(R*X)/(R*X) 

XJR(2)=CSIN(R*X)/( (R*X)**2)-CCOS(R*X)/(R*X) 

XH ( 1 ) =CEXP ( ALPH ) /ALPH 

XH(2)=( (-CEXP(ALPH) )/X)-( (0. , 1 . ) *CEXP{ ALPH) /(X**2 ) ) 
DJ(1)=  COS(X)/X-  SIN(X)/(X**2) 
DJR(l)=CCOS(R*X)/(R*X)-CSIN(R*X)/( (R*X)**2) 
DH(l)=CEXP(ALPH)/X+( (0. ,1. ) *CEXP( ALPH)/(X**2 ) ) 
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QEXT=0 . 0 
QSCA=0 . 0 
C 

C  LI  =  ITERATIONS  FOR  SPHERICAL  FUNCTIONS,  L1=1.3*X 

L1=1.3*X 

IF(Ll.LE.lO)  Ll=10 
DO  1  L=1,L1 
CO=(FLOAT(L) )/(2*L+l) 

COE=(FLOAT(L+l) )/(2*L+l) 

COEF= ( FLOAT ( 2  *L+ 1 ) ) /X 
COEFR=(FLOAT(2*L+l) )/(R*X) 

XJ(L+2)=COEF*XJ(L+l)-XJ(L) 

XJR ( L+2 ) =COEFR*XJR ( L+1 ) -XJR ( L ) 

XH ( L+2 ) =COEF*XH ( L+1 ) -XH ( L ) 

DJ(L+l)=CO*XJ(L)-COE*XJ(L+2) 

DJR ( L+1 ) =CO*XJR ( L ) -COE*X JR ( L+2 ) 

DH ( L+1 ) =CO*XH ( L ) -COE*XH ( L+2 ) 

P(L)=X*XJ(L+1) 

PR(L)=(R*X)*XJR(L+1) 

DP ( L) =XJ ( L+1 ) +X*D J ( L+1 ) 

DPR(L)=XJR(L+1)+(R*X)*DJR(L+1) 

E(L)=X*XH(L+1) 

DE(L)=XH(L+1)+X*DH(L+1) 

A(L)=( (DPR(L)*P(L) )-(R*PR(L)*DP(L) ) )/( (DPR(L)*E(L) )-(R*PR(L)*DE(L 

B(L)  =  ( (R*DPR(L)*P(L) ) - (PR(L) *DP( L) ) )/( (R*DPR( L) *E ( L) ) - (PR( L) *DE  (L 
1))) 

C(L)=A(L)+B(L) 

GAMMA=CABS ( A ( L ) ) 

BETA=CABS(B(L) ) 

XK(L)  =  (FLOAT(2*L+l)  )*(2./(X**2)  )  *(  (GAI'IMA**2)  +  (BETA**2  )  ) 
Xy(L)=(FLOAT(2*L+l) )*(2./(X**2))*(REAL(C(L) ) ) 

QEXT=QEXT+XY(L) 

QSCA=QSCA+XK(L) 

1  CONTINUE 

QABS=QEXT-QSCA 
DQ(M,1)=D 
DQ(M,2)=QEXT 

PRINT  22,X,DQ(M,1) ,QABS,DQ(M,2) 

22  F0RMAT(1X,F8.3 ,6X,F7.3,9X,F6.3,9X,F6.3) 

D  AND  QEXT  FOR  MONODISPERSE  PUT  INTO  PLOTTING  ARRAY 
EVERY  FOURTH  POINT  CALCULATED  IS  PUT  INTO  ARRAY-  ALLQ 

XM=M 

AREL=XM/4 . 0 
MINT=M/4 
XMINT=MINT 

IF  (AREL.NE.XMINT)  GO  TO  60 
IF  (D.GT.1.20)  GO  TO  60 
C  NUMBER  OF  POINTS  FOR  PLOTS  1  AND  3  IS  NPTSl 

NPTS1=MINT 
ALLQ ( MINT+1 , 1 , NI ) =D 
ALLQ(MINT+1 , 3 ,NI)=QEXT 
ALLQ ( MINT+1 , 1 , NI ) =D 
ALLQ (MINT+1 , 3 , NI ) =QEXT 
60  CONTINUE 

C  INCREMENT  SIZE  PARAMETERS  D=DIAMETER 

D=D+DELD 
DI=DI+DELDI 
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M=M4-1 

2  CONTINUE 

C  PRINT  36 

36  FORMAT  ( IH  ) 

DELD=2.0*DELD 
DELDI=2.0*DELDI 

3  CONTINUE 
PRINT  106,NDATA 
PRINT  235 

235  FORMAT('  MIE  COEFFICIENTS  FOR  LOG-NORMAL  DUST  DISTRIBUTION') 

PRINT  236, SIGMA 

236  FORMAT  ('  SIGMA=' , F4 . 2 ) 

PRINT  237,RX 

237  FORMAT('  RELATIVE  REFRACTIVE  INDEX=' , 2F7 . 4 ) 

PRINT  30, WAVE 

CALL  QAVG ( M , NI , SIGMA , NPTS2 , SD3 2 , SDL3  2 , ND ) 

4  CONTINUE 
PRINT  106,NDATA 

7  CONTINUE 
PRINT  500 

500  FORMAT {'  QBAR  RATIOS  FOR  3  WAVELENGTHS') 

PRINT  236, SIGMA 
PRINT  237, RX 

PRINT  509,WAV(1) ,T1,WAV(2) ,T2,WAV(3) ,T3 

509  FORMAT{'  EXP.  TRANSMITTANCES : ' , 3 ( 5X , 'T( ' , F6 . 4 , ' )=' , F6 . 4 ) ) 

PRINT  510,WAV(3) ,WAV(1) ,WAV(3) ,WAV(2) ,WAV(2) ,WAV(1) 

510  FORMAT(/2X,'  D32  ' , 3 ( 5X , 'QBAR( ' ,F6 . 4 , ' ) /QBAR( ' , F6 . 4 , ' ) ' ) , 5X , ' SUM 
+OF  SQUARED  DEVIATIONS'/) 

K=ND+1 

DO  600  N=2,K 

RAT1(N)=ALLQ(N,4,3)/ALLQ(N,4,1) 

RAT2(N)=ALLQ(N,4,3)/ALLQ(N,4,2) 

RAT3(N)=ALLQ(N,4,2)/ALLQ(N,4,1) 

RAT4(N)=(RAT1(N)-L0G(T3)/L0G(T1) )**2 
RAT5(N)=(RAT2(N)-LOG(T3)/LOG(T2) )**2 
RAT6(N)=(RAT3(N)-LOG(T2)/LOG(Tl) )**2 
SSD(N)=RAT4(N)+RAT5(N)+RAT6(N) 

PRINT  650,DIAM32(N) ,RAT1(N) ,RAT2(N) ,RAT3(N) ,SSD(N) 

600  CONTINUE 

650  FORMAT(lX,F6,3,13X,F7.4,23X,F7.4,23X,F7.4,23X,F9.4) 

PRINT  350,NDAT 

350  FORMAT  (/'**************************  END  OF  DATA  SET  NUMBER', 13/) 
SMALLEST=2 
DO  700  N=2,K 

IF  (SSD(N) .LT.SSD(SMALLEST) )  SMALLEST=N 

700  CONTINUE 

PRINT  701, RDUST, SIGMA, DIAM32(SMALLEST) , SSD ( SMALLEST) 

701  FORMAT  ( IX , 2F7 . 4 , ' I ' , 5X , F4 . 2 , 6X , F5 . 2 , 6X , F13 . 8 ) 

NDAT=NDAT+1 
RDUST=RDUST+STEPR 

8  CONTINUE 
RDUST=RDUST1 
SIGMA=SIGMA+STEPS 

9  CONTINUE 
END 


SUBROUTINE  QAVG ( M , NI , SIGMA , NPTS 2 , SD3  2 , SDL3 2 , ND ) 
BY  KENNETH  CASHDOLLAR,  1976 
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C  QAVG  CALCULATES  D32  AND  QBAR  FOR  A  GIVEN  QEXT  VS  D  AND 

C  A  GIVEN  LOG-NORMAL  SIZE  DISTRIBUTION 

C 
C 

C  INPUT  DATA  =  ND, SIGMA, D32,DEL32 

C 

C 

COMMON  ALLQ(100,4,3) 

COMMON/ADQ/DQ (500,2) 

COMMON  DIAM32(100) 

C  DQ(I,1)=D,  DQ(I,2)=QEXT 

C  NI=WAVELENGTH  OR  REFRACTIVE  INDEX  CHANGE 

C  ND  IS  NUMBER  OF  D32,QBAR  CALCULATIONS 

C  SIGMA  IS  THE  GEOMETRIC  STANDARD  DEVIATION 

C  D32  IS  THE  SURFACE  WEIGHTED  MEAN  DIAMETER 

C  DEL32  IS  INCREMENT  FOR  D32 

C  DIST  IS  SURFACE  WEIGHTED  LOG-NORMAL  SIZE  DISTRIBUTION 

C  M  IS  THE  MAXIMUM  NUMBER  OF  POINTS  IN  Q-ARRAY  FOR  MONODISPERSE 

C 

C  PRINT  240 

C  240  FORMAT(3X, 'D32' ,8X, 'DMIN' ,7X, 'DMAX' ,5X, 'NMBR' ,6X, 'DlOX' ,7X, 'D32X' 
C  1,7X,'QBARV) 

C 

D32=SD32 
DEL32=SDL32 

NUMBER  OF  POINTS  FOR  PLOTS  2  AND  4  IS  NPTS2 
NPTS2=ND+1 
DO  17  NN=1,ND 
DMIN=D32/SIGMA**3 . 

DMAX=D 3 2 * S I GMA* * 2 . 5 
AS=ALOG ( SIGMA ) *  *  2 . 

DG=EXP ( ALOG (D32)-2.5*AS) 

QSUM=0 . 0 
DSUM=0 . 0 
XSUM=0 . 0 
DHOLD=0 . 0 
DNSUM=0 . 0 
XNSUM=0 . 0 
NMBR=0 
MM=M-2 

DO  13  N=1,MM 

IF  (DQ(N,1) .LT.DMIN)  GO  TO  13 
IF  (DQ(N,1) .GT.DMAX)  GO  TO  15 
IF  (NMBR.GT.O)  GO  TO  242 
DMINZ=DQ(N,1) 

242  CONTINUE 

DD=DQ(N+1,1)-DQ(N,1) 

AA=(ALOG(DQ(N+l,l)/DG) )**2/(2.0*AS) 
DIST=EXP(-AA)*DQ(N+1,1)/SQRT(6.2832*AS) 

QSUM=QSUM+ ( DIST*DQ ( N+1 , 2 ) +DHOLD*DQ ( N , 2 ) ) *DD/2 . 0 
DSUM=DSUM+ ( DIST*DQ ( N+1 , l)+DHOLD*DQ( N , 1 ) ) *DD/2 . 0 
XSUM=XSUM+ ( DIST+DHOLD ) *DD/2 . 0 

XNSUM=XNSUM+ ( DI ST/DQ ( N+1 , 1 ) *  *  2+DHOLD/DQ (N,1)**2)/2.0*DD 
DNSUM=DNSUM+ ( DIST/DQ ( N+1 , 1 ) +DHOLD/DQ ( N , 1 ) ) /2 . 0*DD 
XHOLD=DQ(N+l,l) 

DHOLD=DIST 
NMBR=NMBR+1 
13  CONTINUE 
15  CONTINUE 
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QBAR=QSUM/XSUM 
D32X=DSUM/XSUM 
D10X=DNSUM/XNSUM 
ALLQ ( NN+1 , 2 , NI ) =D3  2X 
ALLQ ( NN+1 , 4 , NI ) =QBAR 
DIAM32(NN+1)=D32 

PRINT  241,032, DMIN Z , XHOLD , NMBR ,D10X,D32X, QB AR 
241  F0RMAT(1X,3(F6.3,5X) , 13 , 5X , 2 ( F6 . 3 , 5X) , F7 , 4 , 5X) 
D32=D32+DEL32 
17  CONTINUE 
RETURN 
END 
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APPENDIX  F.  LINEAR  INTERPOLATION  IN  LOTUS  1-2-3  RELEASE  4 


A  formula  was  written  to  perform  linear  interpolation  in  Lotus®  1-2-3®  Release  4 

for  Windows'™.  The  formula  used  the  diameter  value  in  column  B  (cell  address  BIO 

below)  to  choose  the  correct  extinction  coefficient  Q  from  a  separate  1-2-3  file  named 

"qext.wk4".  The  file  "qext.wk4"  had  a  named  range  QEXT  containing  seven  columns  of 

data:  diameter  and  six  columns  of  Q's  for  different  wavelengths  and  indices  of  refraction. 

The  formula  below  contains  the  offset  for  the  first  column  of  Q's. 

@  INDEX(«qext.  wk4»$QEXT,  1 ,  @MATCH(B  1 0,«qext.  wk4»$QEXT,  1 ))+(( 
@INDEX(«qext.wk4»$QEXT,  1 ,  @MATCH(B  10,«qext.wk4»$QEXT,  1)+1)- 
@  INDEX(«qext.  wk4»$QEXT,  1 ,  @MATCH(B  1 0,«qext.  wk4»$QEXT,  1 )))/( @ 
INDEX(«qext.wk4»$QEXT,0.@MATCH(B10,«qext.wk4»$QEXT,l)+l)-@I 
NDEX(«qext.wk4»$QEXT,0,@MATCH(B10,«qext.wk4»$QEXT,l)))*(B10 
-  @INDEX(«qext.wk4»$QEXT,0,  @MATCH(B  10,«qext.wk4»$QEXT,  1)))) 

Linear  interpolation  can  best  be  described  by  an  example.  Suppose  di  and  d2  are 
increasing  diameters  with  corresponding  extinction  coefficients  Qi  and  Q^.  Then  the  Q 
for  a  diameter  d  between  d^  and  d^  is: 

Q  =  Qx  +  {^^^*[d-d,]^  (F.l) 


The  reader  is  directed  to  the  User's  Guide  for  further  information  and  the  ©function  help 
in  the  software  for  an  explanation  of  ©INDEX  and  ©MATCH. 


83 


84 


APPENDIX  G.  DATA  PRINTOUTS,  PLOTS,  AND  TABLES 
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Figure  10.  Malvern  Data  from  PS-1  Firing  3/7/95. 
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Figure  11.  PDPA  Data  from  PS-1  Firing  3/2/95. 
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Figure  13.  P  DP  A  Data  from  PS- 1  Firing  3/4/95. 
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Figure  14.  LABTECH  Notebook  Data  Acquisition  Program  Plot  of  PS-1  Firing  2/25/95 
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Figure  15.  Normal  Plot  of  PS-1  Firing  2/25/95. 
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Figure  16.  Bare  Light  (No  Probe)  PS-1  Firing  3/5/95. 
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Figure  19.  Scanning  Electron  Microscope  Pictures  of  AljOj  Collected  in  the  Probe. 
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